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ABSTRACT 


The  effect  of  refraction  on  radar  and  ballistic  camera  tracking  accuracy 
is  presented  against  a  background  of  operational  requirements  and  limitations. 

A  method  to  compute  refraction- induced  errors  using  closed  form  expressions 
is  presented  together  with  a  comparison  of  these  errors  against  ray  tracing 
results.  For  a  hypothetical  missile  trajectory  it  is  shown  that  the  closed  form 
solutions  using  an  exponential  model  atmosphere  are  in  good  agreement  with 
the  ray  tracing  results  using  radiosonde  data.  The  analysis  also  shows  that 
the  closed  form  expressions  can  be  used  to  correct  ballistic  camera  tracking 
errors.  In  this  case  a  modified  form  of  the  exponential  model  is  used  to  de¬ 
scribe  the  optical  refraction  profile.  Variations  in  refractivity  in  time  and 
space  are  examined  for  their  effect  on  tracking  accuracy.  The  analysis  shows 
that  time  variations  of  the  refractivity  profile  near  the  radar  can  produce 
significant  variations  in  the  elevation  angle  error.  Spatial  variations  far 
from  the  radar  are  restricted  in  their  effect  and  because  of  earth  curvature 
only  variations  near  the  radar  are  significant.  For  tracking  below  five  de¬ 
grees  elevation  angle  the  tracking  errors  can  best  be  determined  from  ray 
tracing  analysis  with  real-time,  corrected,  radiosonde  data.  However,  the 
ability  to  track  much  below  five  degrees  is  shown  to  be  impaired  by  multipath 
propagation  effects.  Recommendations  are  given  for  real-time  evaluations 
of  the  ray  tracing  and  closed  form  calculations  of  tracking  errors.  The 
advantages  of  the  simplicity  of  application  of  the  closed  form  equations  for 
real-time  corrections  is  stressed  in  the  report. 
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INTRODUCTION 


A  refraction  investigation  was  carried  out  at  the  Air  Force  Western 
Test  Range  to  determine  the  tracking  errors  which  would  be  produced  by  the 
troposphere.  Refraction  measurements  were  made  using  a  Crain  microwave 
refractometeri<smounted,  together  with  other  instrumentation,  in  a  USAF 
C131B  Convair.  This  facility  is  described  in  detail  together  with  the  ground 
based  analysis  programs  in  the  Appendix. 

Additional  refraction  data  were  obtained  from  rawinsonde  measure¬ 
ments  launched  from  North  Vandenberg,  the  Boathouse,  San  Nicolas  Island, 
and  Point  Mugu.  The  Pacific  Missile  Range  (PMR)  provided  radar  tracking 
and  meteorological  support  on  many  occasions  when  flights  were  conducted 
in  their  range. 

The  initial  effort  began  in  January  1967  through  June  1967,  at  which 
time  the  contract  was  with  the  Mitre  Corporation,  Bedford,  Massachusetts. 

During  this  period  an  analysis  of  the  effects  of  refraction  errors  on  bilatera- 
tion  and  trilateration  tracking  errors  was  presented  to  AFWTR  in  April  1967. 3 

In  June  1968,  the  refraction  investigations  were  reinitiated  at  AFWTR 
under  the  direction  of  the  Syracuse  University  Research  Corporation.  The 
aircraft  and  meteorological  support  facilities  were  again  made  available  on 
this  latter  effort. 

This  report  presents  an  analysis  of  the  refraction  conditions  affecting 
the  radar  tracking  systems  which  support  the  AFWTR  firing  program.  These 
trackers  are  located  at  Tranquillan  Peak,  San  Nicolas  Island,  and  Point  Mugu. 

The  fourth  station  at  Point  Pillar  (near  San  Francisco)  is  not  included  in  the 
investigation  described  herein.  It  was  intended  to  complete  the  program  in 
the  AFWTR  sector  at  Vandenberg  before  initiating  refraction  studies  at 
Point  Pillar. 

In  an  investigation  of  this  kind,  very  large  quantities  of  data  are  col¬ 
lected,  not  only  from  the  aircraft  measurements  but  also  from  the  meteorological 
support  facilities  at  AFWTR.  This  total  compendium  of  data  are  carefully  des¬ 
cribed  and  stored  at  Syracuse  in  the  event  that  further  detailed  analysis  is  re¬ 
quired.  It  is  not  intended  that  the  report  describe  this  total  mass  of  data  but 
rather  selected  situations  are  presented  which  illustrate  the  nature  of  the 
problem  and  from  which  some  definite  conclusions  and  recommendations  can 
be  made. 
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A  concerted  effort  is  made  to  relate  this  investigation  to  previous  re¬ 
fraction  studies  carried  out  at  PMR  and  AFWTR.  In  this  way,  the  present 
investigation  and  analysis  represents  an  extension  of  these  former  programs 
and  takes  advantage  of  the  extensive  ray  tracing  comparisons  reported  earlier 
by  Gardner  .4 

Several  methods  are  discussed  which  can  be  used  to  correct  range  and 
elevation  angle  tracking  errors  together  with  their  relative  advantages.  From 
a  description  of  the  spatial  and  temporal  variations  of  the  radio  refractivity 
the  magnitudes  of  residual  errors  are  calculated  and  their  effect  in  the  system 
is  related  to  operational  procedures 

Since  ballistic  tracking  cameras  experience  tracking  errors  due  to  the 
dry  term  of  refractivity,  a  brief  discussion  of  these  errors  is  presented. 
Consideration  is  given  to  the  use  of  an  exponential  function  to  correct  these 
tracking  errors. 

A  method  to  correct  velocity  measurement  errors  is  presented.  This 
investigation  was  not  intended  to  be  the  major  effort  during  this  contract 
period;  therefore,  the  presentation  is  brief  but  demonstrates  a  correction 
method  which  could  have  considerable  potential. 

The  most  significant  part  of  the  presentation  herein  concerns  the  cal¬ 
culation  of  tracking  errors  from  simple,  closed,  functions.  This  correction 
method  requires  that  the  vertical  variation  of  refractivity  be  represented  by 
an  exponential  model.  In  comparison  with  Gardner's  analysis  it  is  shown 
that  this  simple  method  can  be  used  to  correct  tracking  errors  for  a  wide 
range  of  conditions 

Finally,  recommendations  are  made  for  evaluating  the  error  correc¬ 
tion  programs,  initially  using  post-flight  data  to  be  followed  by  the  application 
of  selected  correction  techniques  during  live  firings.  Methods  are  also  pre¬ 
sented  to  acquire  and  apply  appropriate  meteorological  data  in  the  correction 
programs 
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SECTION  I 


THE  EFFECT  OF  TROPOSPHERIC  REFRACTION  ON  ELECTRO¬ 
MAGNETIC  PROPAGATION 


Before  discussing  refraction-induced  errors  and  correction  techniques 
in  detail,  it  is  useful  to  review  the  nature  of  these  errors  and  general  methods 
for  correction  which  presently  exist. 


The  index  of  refraction  in  the  troposphere  is  given  by6 


=  r77 .  6  4810  e  vl 

n  =  1  +  10  6  [— =—  (P  +  - = - )] 


(1) 


where  P  =  the  pressure  of  a  volume  of  air  (mb) 

T  =  the  temperature  of  the  air  (°  K) 
e  =  the  water  vapor  pressure  in  the  volume  (mb) 


The  variable  part  of  the  index  is  called  the  refractivity  and  is  designated 


nt  = 


77. 6 
T 


<F  +  «^) 


(2) 


where  the  dry  term  is 


ndry 


77  6  P 
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NWET  " 


77.  6 
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4810  e 

V  rj-  / 


(3) 

(4) 


The  velocity  of  propagation  of  electromagnetic  energy  is  given  by 


v  =  c  /  n 


(5) 


where  c  is  the  velocity  of  light  or  the  velocity  of  electromagnetic  propaga¬ 
tion  in  a  vacuum 


From  propagation  analysis  it  can  be  shown  that  at  a  point  in  space  the 
reciprocal  of  the  radius  of  curvature  of  an  electromagnetic  signal  (ray  theory) 
is  a  function  of  the  gradient  of  the  logarithm  of  refraction,  n  6  The  most 
significant  and  sustained  gradients  (variation  of  n  with  distance)  occur  in  the 
vertical  direction  and  the  curvature  causes  the  signal  to  propagate  towards 
the  region  of  increasing  index 

For  radio  frequencies  the  index  of  refractivity  which  affects  the  velocity 
and  direction  of  propagation  is  given  by  the  total  term  (Equation  (2))s  With 
the  exception  of  water  vapor  resonance  at  22  GHz  and  oxygen  resonance  at 
60  GHz  the  index  is  essentially  frequency  independent  up  to  72  GHz  ,8>7 
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In  the  visible  optical  region  the  refractivity  is  expressed  by  only  the 
dry  term  (Equation  (3))  and  the  effect  of  water  vapor  can  be  neglected  6 

In  both  the  radio  and  optical  frequency  range  the  magnitude  of  re¬ 
fractivity  tends  to  decrease  exponentially  with  height  due  to  the  overwhelming 
influence  of  the  exponentially  height-dependent  pressure  term  8  However,  at 
radio  frequencies  the  variability  of  the  water  vapor  pressure  near  the  earth’s 
surface  (generally  under  25  000  feet)  can  cause  large  variations  from  an 
average  exponential  behavior 

Direct  measurements  of  refractivity  can  be  made  with  microwave  re- 
fractometers  and  pressure,  temperature,  and  water  vapor  instruments.  From 
these  measurements  a  spatial  pattern  of  the  refractivity  can  be  developed 
with  the  most  sustained  variations  occunng  in  the  vertical  direction.  In  most 
cases  the  refractivity  structure  is  defined  in  the  vertical  dir ection,  and  it  is 
assumed  that  the  measurements  can  be  applied  at  other  locations  over  the 
earth's  surface  thereby  assuming  a  spherically  symmetrical  atmosphere. 

Figure  1  defines  the  geometry  associated  with  ray  tracing  analysis  of 
the  effect  of  the  vertical  variation  of  refractivity  on  angle  and  range  errors. 
Using  a  ray  theory  approach  the  signal  is  transmitted  to  and/or  received 
from  a  target  at  T.  The  apparent  elevation  angle  of  the  object  is  90,  whereas 
the  true  elevation  angle  is  fl0.  Due  to  the  downward  curvature  of  the  ray  an 
elevation  angle  error,  e,  is  produced.  In  theory,  if  the  refractivity  of  the 
medium  was  accurately  known  the  magnitude  of  the  elevation  angle,  c,  could 
be  determined  by  ray  tracing  analysis  for  given  geometrical  coordinates  at 
each  end  of  the  ray  (Ray  tracing  methods  wrill  be  discussed  later.  ) 

The  true  geometrical  range  to  the  target  is  R0.  The  measured  range 
to  the  target  is  along  the  electrical  path,  Re .  By  measuring  the  time  for  a 
signal  to  travel  between  the  end  points  of  the  ray  (echo  ranging)  and  assuming 
that  the  velocity  of  propagation  is  c  (velocity  in  vacuum)  two  errors  result. 

The  first  error  results  from  the  fact  that  the  geometrical  distance  along  Re, 
calling  this  distance  Rg,  is  longer  than  the  distance  along  R0.  The  target, 
therefore,  appears  to  be  farther  away  than  it  actually  is.  The  second  error 
results  from  the  fact  that  the  velocity  of  propagation  is  less  than  the  velocity 
in  vacuum  by  an  amount  n"1  (Equation  (5))  at  every  point  along  the  ray.  The 
delay  produced  by  this  reduced  propagation  velocity  along  the  ray  causes  a 
retardation  error  in  the  range  measurement  which  will  be  designated,  AR . 

The  distance  along  Re,  measured  by  the  round-trip  travel  time  of  a  signal, 
is  therefore  too  large  due  to  velocity  retardation.  Ray  tracing  analysis  has 
shown  that  the  geometrical  or  first  error,  Rg  -  R0,  is  a  second  order  error 
compared  to  the  retardation  error,  AR/  Therefore,  when  range  error  is 
mentioned  herein  we  are  refering  to  the  retardation  error,  AR. 
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FIGURE  1.  RAY  GEOMETRY 
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SECTION  II 


THE  DETERMINATION  OF  TFACKING  ERRORS  BY  RAY¬ 
TRACING  ANALYSIS 

The  determination  of  range  and  directional  errors  can  be  carried  out 
with  great  accuracy  using  digital  computer  facilities.  The  usual  method  is 
to  assume  that  a  given  refractivity  profile  (r efractivity  versus  height)  repre¬ 
sents  a  spherically  symmetrical  condition  over  the  earth's  surface.9  The 
troposphere  is  broken  up  into  small  enough  height  sections  such  that  the 
index  gradient,  dn/dh,  can  be  consider  ed  to  be  constant  within  any  given 
section.10  Since  n  will  vary  with  height  the  magnitudes  of  these  gradients  will 
generally  differ  for  each  height  section  From  optical  ray  theory,  the 
incremental  bending  of  a  ray  can  be  calculated  as  it  passes  through  a  particular 
section.6  Similarity  the  retardation  produced  in  the  section  can  be  determined 
By  summing  these  effects  it  is  possible  to  calculate  the  total  bending,  r,  im¬ 
pressed  on  the  ray  over  the  path  Re  and  to  determine  the  retardation  error, 

AR,  Once  the  path  of  the  ray  has  been  determined  in  this  manner  all  other 
parameters  are  available  from  these  calculations,  such  as  the  true  range,  Rq, 
the  elevation  angle  error,  e,  and  the  doppler  tracking  error  angle,  6  (refer¬ 
ence  Figure  1).  Discussions  on  the  effect  of  this  angle,  6,  on  velocity 
measurement  errors  will  be  presented  later  , 

A  comparison  of  ray-tracing  programs  used  within  various  ranges 
shows  that  in  terms  of  mathematical  precision  there  is  essentially  no  dis¬ 
agreement,  The  7030  program  designated  Refchex  used  on  th*»  s  contract 
has  been  checked  against  Gardner's  REF  PAG  and  Bean  and  Thayer's  program,  1? 
showing  that  there  are  no  significant  differences.  Therefore,  the  accuracy 
available  with  present  digital  ray-tracing  programs  is  not  a  limitation  to  the 
determination  of  tracking  errors  for  a  prescribed  index  profile. 

There  are,  however,  certain  considerations  which  affect  the  useful¬ 
ness  of  ray-tracing  techniques  for  the  real-time  correction  of  tracking  data 

(i)  In  the  first  place,  can  a  refractivity  profile  be  determined  and 
used  in  the  error  correction  program  in  a  short  enough  period 
of  time  such  that  current  conditions  are  represented? 

(ii)  Second,  what  errors  are  introduced  by  the  assumption  that 
this  refractivity  measurement  is  spatially  invariant? 

(m)  Third,  if  the  tracking  errors  are  calculated  by  ray  tracing  with 
this  particular  profile,  how  much  calculation  time  and  com¬ 
puter  storage  is  required  to  cover  the  limits  of  tracking  angles 
and  ranges? 

(iv)  Fourth,  can  this  error  data  be  acquired  from  storage  quickly 
enough  to  correct  real-time  tracking  data? 
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SECTION  III 


THE  SPATIAL  AND  TEMPORAL  VARIATIONS  OF  THE 
REFRACTIVE  INDEX 

Aircraft  measurements  of  the  refractive  index  profiles  were  made 
from  near  sea  level  to  about  15,  000  feet  and  at  intervals  in  distance  out  to 
200  nautical  miles  from  Vandenberg  AFB  and  Point  Mugu.  Figure  2  shows 
a  map  of  the  area  with  the  tracks  generally  flown  along  a  westward  outbound 
bearing.  As  shown,  spirals  were  repeated  on  the  inbound  flights  at  the  same 
coordinates  as  outbound  to  determine  the  changes  which  occurred  in  time. 

A  rawinsonde  launch  was  scheduled  to  coincide  with  the  first  aircraft 
spiral  on  the  outbound  leg.  Although  the  aircraft  instrumentation  was  cali¬ 
brated  against  psychometric  measurements  on  the  ground,  prior  to  the  flight 
the  rawinsonde  profile  provided  a  check  on  the  aircraft  measurements  in  flight. 

1 .  Aircraft  and  Rawinsonde  Measurements 

Figures  3a,  b,  and  c  show  typical  comparisons  of  aircraft  and  rawin¬ 
sonde  measurements  As  the  aircraft  ascends  it  circles  about  a  vertical  axis 
which  is  generally  three  miles  diameter  at  altitudes  below  3,  000  feet  and  about 
four  miles  diameter  at  15,000  feet.  For  this  reason,  much  of  the  fine 
structure  shown  on  the  aircraft  profile  can  be  attributed  to  horizontal  varia¬ 
tions  It  is  pertinent  to  later  discussions  to  note  on  Figure  3a  that  the  rawin¬ 
sonde  does  not  respond  to  an  inversion  which  was  measured  by  the  aircraft. 

The  aircraft  was  about  10  miles  seaward  from  the  rawinsonde  launch  site. 

In  this  connection,  a  difficulty  is  presented  by  the  fact  that  the  winds 
tend  to  be  westwardly  in  which  case  rawinsondes  launched  from  the  coastal 
sites  are  blown  inland.  Figure  4  shows  representative  trajectories  of  rawin¬ 
sondes  launched  from  North  Vandenberg  and  the  Boathouse.  Figure  5  shows 
a  plot  of  wind  direction  versus  altitude  from  eight  rawinsonde  launches,  four 
from  North  Vandenberg  and  four  from  the  Boathouse.  In  this  instance,  the 
wind  has  shifted  to  the  northwest  over  the  period  29-31  May  1967  In  general, 
meteorological  measurements  provided  by  AFWTR  and  PMR  over  the  test 
periods  show  the  wind  prevailing  from  west  to  northwest.13  It  was  not  possible 
for  the  aircraft  to  track  the  rawinsonde  since  visibility  near  the  launch  areas 
was  usually  obstructed  by  fog  and/or  haze. 

Due  to  the  drift  on  most  rawinsonde  launches  their  measurements  are 
not  in  the  undisturbed  maritime  air  mass  which  is  west  of  the  Tranquillan 
Peak  and  Point  Mugu  tracking  stations  Comparisons  with  aircraft  measure¬ 
ments  have  not  shown  any  significant  local  differences  except  in  situations 
when  a  strong  maritime  inversion  is  present.  Due  to  the  inland  drift  of  the 
rawinsondes,  the  intensity  of  these  inversions  can  frequently  be  underestimated. 
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FIGURE  2.  FLIGHT  TEST  AREA 
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FIGURE  3a.  A  COMPARISON  OF  AIRCRAFT  AND  RAWINSONDE 
REFRACTIVITY  PROFILES 
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FIGURE  3c.  A  COMPARISON  OF  AIRCRAFT  AND  RADIOSONDE 
REFRACTIVITY  PROFILES 
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FIGURE  5.  BOATHOUSE  AND  VANDENBERG  RADIOSONDES 


San  Nicolas  Island  launches  were  expected  to  represent  fairly  undisturbed 
maritime  conditions.  It  was  initially  anticipated  that  these  launches  could  pro¬ 
vide  a  meaningful  description  of  refractivity  for  use  by  all  the  tracking  stations 
in  the  AFWTR  and  PMR  areas 

Figures  6,  7,  and  8  show  a  comparison  of  profiles  for  simultaneous 
launches.  In  general  there  is  not  good  agreement  until  heights  near  twenty 
thousand  feet  have  been  reached  In  Figure  7,  even  in  the  presence  of  a  strong 
inversion  (vertical  stability)  at  San  Nicolas  and  Point  Mugu  and  a  weak  in¬ 
version  at  Vandenberg  ,there  are  significant  variations  occurred  above  the  in¬ 
version  layert 

Aircraft  measurements  of  the  inversion  characteristics  and  observa¬ 
tions  of  cloud  structure  indicate  there  are  large  spatial  variations  in  water 
vapor  and  visible  water,  respectively.  At  heights  below  ten  thousand  feet 
there  are  large  differences  in  the  magnitudes  of  water  vapor  pressure  between 
stations  Therefore,  evidence  indicates  that  the  San  Nicolas  rawinsonde  can¬ 
not  be  used  to  determine  a  profile  for  general  application  to  all  the  tracking 
stations . 

2.  Aircraft  Measurements  of  the  Radio  Refractivity  in  the  Seaward 

Direction 

From  aircraft  measurements,  it  was  possible  to  get  in  indication  of 
the  spatial  and  temporal  variations  of  radio  refractivity  along  the  seaward  ex¬ 
tension  of  the  range  from  Vandenberg  Figures  9  and  10  show  refractivity 
profiles  obtained  along  a  flight  path  bearing  two  hundred  and  sixty  degrees 
from  Vandenberg.  Profiles  A1#  B1,  and  C1  represent  vertical  spirals  out¬ 
bound  and  B2,  Az  co-located  spirals  inbound  Soundings  made  during  the  climb 
from  the  bottom  of  one  spiral  to  the  top  of  the  next  vertical  spiral  are  desig¬ 
nated  B1-C1  and  B2-Cl  (Figure  9)-  Both  sets  of  data  show  large  variations 
were  measured  between  profiles. 

Comparing  spirals  A1  and  A2  which  are  at  the  same  point  but  different 
in  time,  it  can  be  seen  that  the  characteristics  of  the  inversion  are  changing 
significantly  over  a  period  of  about  three  hours 

Figure  11  derived  from  Figure  9,  shows  the  average  profile  and  the 
variations  which  occur  red  as  a  result  of  both  spatial  and  temporal  changes. 
Figures  12  and  13  show  the  isopleths  of  constant  refractivity  derived  from 
the  soundings  of  Figure  9  Superimposed  on  these  latter  figures  are  curves 
showing  the  radar  horizon  for  elevation  angles  of  zero,  one,  and  two  degrees. 

In  general,  tracking  below  two  degrees  subjects  the  radar  to  sea  reflections. 
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FIGURE  6.  RADIOSONDES  OCTOBER  23,  1967  1715Z 


15 


HEIGHT  (KFT) 


CD 


FIGURE  7.  RADIOSONDES  OCTOBER  23,  1967  1100Z 
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FIGURE  8.  RADIOSONDE  JUNE  13,  1967 
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FIGURE  9.  RADIO  REFRACTIVITY  PROFILES  WEST  FROM 
VANDENBERG,  JUNE  26,  1968 
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Figure  14  shows  a  typical  range-height  plot  from  the  Tranquillan  Peak  radar. 
The  aircraft  made  an  initial  upward  spiral  at  a  point  10  nautical  miles  from 
the  radar  and  terminated  the  spiral  at  15,  000  feet  above  sea  level.  The  air¬ 
craft  then  descended  to  reach  a  coordinate  approximately  85  miles  from  the 
radar  where  a  second  spiral  was  initiated.  Since  no  earth  curvature  or  re¬ 
fraction  correction  was  used,  the  plot  shows  the  top  of  the  second  spiral  to 
be  approximately  10,  000  feet  above  sea  level  whereas  the  true  height  is 
15,  000  feet. 

Using  a  simple  4/3  effective  earth  radius^  the  effect  of  sea  reflections 
on  the  radar  elevation  angle  data  was  pronounced  for  apparent  elevation 
angles  below  1.  5  degrees. 

Figure  15  shows  the  ray  path  geometry  where  the  true  earth  radius 
has  been  replaced  by  an  effective  earth  radius.  The  Tranquillan  Peak  radar 
is  about  2500  feet  above  sea  level.  From  Figure  15  it  is  apparent  that  the 
angle  between  the  direct  and  reflected  signals,  O',  becomes  less  as  the  height 
of  the  antenna  above  the  reflecting  sea  surface  decreases.  For  tracking 
stations  at  Point  Mugu  and  San  Nicolas  Island,  which  are  at  lower  heights 
above  sea  level,  the  reflected  signal  arrives  much  closer  to  the  direction  of 
the  main  beam.  If  the  sidelobe  level  is  not  very  far  down  from  the  main 
beam,  the  behavior  of  these  monopulse  tracker  systems  can  be  severely 
disturbed  in  the  elevation  angle  measurement. 

It  is,  therefore,  reasonable  to  suggest  that  tracking  below  1.  5  degrees 
apparent  elevation  angle  cannot  be  carried  out  with  the  existing  radars  in  the 
AFWTR  system.  Only  when  the  target  is  on  the  sea  surface  would  multipath 
propagation  be  prevented  In  fact,  the  radars  can  track  ships  without  serious 
elevation  angle  fluctuations. 

3 .  The  Effect  of  Time  and  Space  Variations  of  the  Refractivity  Profile 

Let  it  be  assumed  that  Profile  A1,  Figure  9,  represents  an  accurate 
determination  of  radio  propagation  conditions  pertaining  to  the  Tranquillan 
Peak  radar  at  1639  hours.  Restricting  tracking  to  two  degrees  elevation 
angle  or  greater  it  is  reasonable  to  assume  that  Profile  A1  represents  spatial 
conditions  out  to  50  nautical  miles  (reference  Figure  12). 
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(AIRCRAFT  AT  15  KFT  )  TRANQUILLAN  PEAK  RADAR 


o 

o 


o 

CD 


o 

ao 


cr  o 
o  r: 


C/) 


Ui 


<  x 
^  a: 

Q 

< 

a: 


o 

CVJ 


o 


♦  ♦  ♦ - f- 

^  ^  ro  w 


»  >  I - I  < - 1  ■»*■  ♦ 

O  CD  00  h-  CD  in  fO  CVJ 


-*  o 
o 


tux)  aaniinv  i3oavi  Q3ansv3w 


X 

H 

< 

W 

H 

C 

X 

p 

H 

o 

P 

X 

H 

s 

a 

i— i 

w 

< 

p 

c 

ccj 

■>* 

r-H 

w 

oi 

D 

0 

i— i 


24 


A  4816 


T 


FIGURE  15.  MULTIPATH  GEOMETRY 


Suppose  that  the  application  of  this  profile  in  the  refraction  correction 
program  does  not  take  place  until  1931  hours.  By  this  time,  Profile  A2 
(Figure  9)  is  now  assumed  to  be  a  correct  description  of  the  radio  propagation 
conditions.  What  errors  occur  as  a  result  of  the  changes  in  refractivity 
between  1639  and  1931  hours? 

Figure  16  shows  the  difference  in  the  elevation  angle  errors,  as 

a  function  of  the  initial  elevation  angle,  Q0,  for  the  Tranquillan  Peak  radar. 

At  the  smallest  tracking  angle  of  two  degrees  a  0  25  milliradian  error  would 
exist  as  a  result  of  the  above  time  variations  in  the  propagation  conditions. 

Figure  17  shows  that  a  2.25  foot  range  error  difference  (A(AR)), 
would  exist  at  two  degrees  elevation  angle.  Ray  tracing  analyses  were  not 
carried  out  above  12,  000  feet  because  radio  refractivity  variations  above 
this  level  were  insignificant  Since  the  aircraft  measurement  accuracy  is 
rechecked  after  each  flight,  the  above  profile  variations  were  determined  to 
be  due  to  temporal  changes  and  not  due  to  instrumentation  errors  in  the  air¬ 
craft. 


With  one  aircraft  it  is  not  possible  to  obtain  time  invariant  spatial 
measurements;  therefore,  the  information  presented  in  Figures  12  and  13 
includes  some  time  variation  effects.  Let  it  be  assumed  that  the  measure¬ 
ment  of  refractivity,  Profile  Ax  (Figure  9),  is  again  put  into  the  tracking 
system  refraction  correction  program.  Approximately  one  hour  later,  a 
second  measurement,  Profile  B1#  is  obtained  at  50  nautical  miles  range  from 
the  tracking  station  A  comparison  of  these  profiles  on  the  elevation  angle 
error  differences,  Ac,  and  the  range  error  difference,  A(AR),  is  shown  in 
Figures  18  and  19,  respectively  Obviously,  from  Figure  12,  the  only  varia¬ 
tions  affecting  tracking  at  two  degrees  occur  above  7,  000  feet  altitude 

At  two  degrees  initial  elevation  angle  the  elevation  angle  error  dif¬ 
ference  is  0.  15  milliradian  and  the  range  error  difference  is  0.42  feet. 

For  tracking  stations  operating  closer  to  sea  level  the  radio  horizon 
at  two  degrees  elevation  angle  provides  even  greater  restriction  on  the  effect 
of  spatial  variations.  As  shown  by  Figure  12,  for  the  Point  Mugu  radar, 
variations  below  7,  000  feet  would  not  be  significant  for  ranges  beyond  35 
nautical  miles. 

It  is  apparent  that  spatial  variations,  which  occur  largely  at  low 
altitudes,  are  prevented  by  the  earth's  curvature  from  having  a  significant 
effect  on  the  radar  tracking  accuracy.  However,  temporal  effects  can  pro¬ 
duce  large  errors  and  should  be  constrained  by  incorporating  refractivity 
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VARIATIONS  OF  REFRACTIVITY 
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FIGURE  17.  RANGE  ERRORS  PRODUCED  BY  TIME  VARIATIONS 
OF  REFRACTIVITY 
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VARIATIONS  OF  REFRACTIVIT  Y 
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VARIATIONS  OF  REFRACTIVITY 


profile  information  into  the  tracking  system  as  close  as  possible  to  the  time 
of  application.  It  is  well  known  that  ray  bending  near  the  radar  produces 
larger  height  errors  than  bending  which  occurs  far  from  the  radar.  This  is 
simply  because  the  height  error  is  the  product  of  the  angle  error  and  the 
distance  from  this  bending  disturbance  to  the  target. 

In  summary,  the  effects  of  temporal  variations  near  the  radar  site 
are  most  effective  in  producing  elevation  angle  errors  and,  therefore,  height 
errors.  For  tracking  above  two  degrees  elevation  angle  spatial  variations 
have  a  secondary  effect  in  the  tracking  accuracy.  Therefore,  evidence  in¬ 
dicates  that  to  determine  refractivity  profiles  for  the  correction  of  refraction- 
induced  errors,  measurements  should  be  made  close  to  the  radar  site  and 
along  a  bearing  of  the  expected  azimuth  tracking  sector.  Finally,  the  re¬ 
fractivity  data  should  be  collected  and  used  as  close  as  possible  to  the  time 
of  tracking  application  in  order  to  constrain  temporal  errors. 
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SECTION  IV 


A  COMPARISON  OF  METHODS  TO  CORRECT  REFRACTION 
ERRORS  AT  RADIO  FREQUENCIES 

At  radio  frequencies  we  are  concerned  with  the  behavior  of  the  total 
index  of  r efr activity ,  N^p,  given  by  Equation  (2).  From  the  analysis  and 
discussion  of  spatial  refractivity  effects  it  will  be  assumed  that  a  profile  can 
be  measured  which  represents  a  description  of  spherically  stratified  radio 
propagation  conditions 

The  first  method  to  correct  refraction  induced  errors  uses  a  ray¬ 
tracing  analysis.  Range  and  angle  errors  are  calculated  for  tracking  eleva¬ 
tion  angles  and  radar  ranges  corresponding  to  the  trajectory  of  a  hypothetical 
missile  launch.  This  analysis  was  carried  out  by  Gardner4  and  his  ray¬ 
tracing  data  will  be  used  to  compare  with  an  alternative  method  to  calculate 
tracking  errors,  desc  rib  ed  below. 

1.  Refraction  Error  Calculations  in  an  Exponential  Atmosphere 

Let  it  be  assumed  for  the  moment  that  the  vertical  decrease  of  re¬ 
fractivity  can  be  expressed  in  terms  of  an  exponential  function8  relative  to 
the  tracking  station,  where 

N  =  Ns  e'Chs  (6) 

and  Ns  =  the  surface  refractivity  at  the  tracking  station 

hs  =  the  height  above  the  elevation  of  the  station  (km)  and 

c  -  the  appropriate  exponential  decay  constant  (km-1). 

From  Figure  1  it  is  apparent  that  the  height,  h,  of  any  point,  P,  on 
the  slant  path,  R0,  can  be  expressed  in  terms  of  the  earth  radius,  rQ,  the 
station  elevation,  Hs,  the  true  elevation  angle,  80,  and  the  geometrical  dis¬ 
tance,  S,  of  the  point,  P,  from  the  station.  In  this  way  the  variation  of  N 
(Equation  (6)),  can  be  expressed  in  terms  of  the  distance,  S,  along  R0  which 
is  now  taken  to  be  the  variable. 

From  ray  theory  it  can  be  shown  that  the  range  error  due  to  retarda¬ 
tion  AR  is  given  by 

Rg 

AR  ~  10"6  J  N.  ds  (7) 

o 

where,  from  previous  discussion  (reference  Section  J)  the  geometrical 
error,  Rg-R0,  is  neglected.7  Since  Rg  and  R0  are  very  nearly  equal,  it  is 
convenient  to  rewrite  Equation  (7)  where 
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(8) 


Ro 

AR  ~  10"6  f  N.ds 
o 

and  integration  takes  place  along  the  slant  path,  R0,  rather  than  along  the 
ray  path,  R^.  Since  N  can  be  expressed  as  a  function  of  distance  along  R0l 
the  integral  (Equation  (8))  can  be  solved  exactly. 

Freeman's  development  of  this  method14  and  solution  showed  that  the 
range  error  due  to  retardation  could  be  expressed  in  terms  of  error  functions. 
Recognition  should  be  given  to  Thayer's  earlier  investigations  which  simulated 
these  further  analyses.15 


Freeman's  solutions  are  given  by  the  following  equations  where 

nr_ 

2r, 


AR  =  — - N-s  -^-|-  e-P---'g'  ~  [erf  (R0  cos  80  /-£—  +  g)  -  erf  (g)] 

c  sin  p0 


where  erf  e  = 


g=vT  J  exp  (-: 
g  -  tan  j30 


dx 


o 

C  Vr 


and  it  is  assumed  that  r0>>Hs. 

For  our  present  discussion,  we  are  dealing  with  tracking  stations 
located  near  sea  level  so  the  above  equations  apply. 


(9) 

(10) 

(H) 


In  a  separately  reported  analysis,  Freeman  showed  that  the  overall 
path  bending,  T,  (Figure  1)  was  given  by  16 

0r>  Nr 

T  ~  1 0 -s  f  cot  0  dN  (12.  1) 

0O»NS 


IQ-6  g  exp  (gs) 

tan  0O 


[erf  (R0  cos 


+  g)  -  erf  (g)] 


(12.2) 


where  9r  and  Nr  are  defined  at  the  target  location. 


These  equations  were  found  to  be  in  error  since  the  values  of  re- 
fractivity  along  Rq  are  larger  than  they  would  be  along  the  higher  path,  Re. 
Therefore,  an  empirically  determined  adjustment  was  made  to  give  agree¬ 
ment  with  published  data  and  for  targets  above  the  troposphere. 
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To  make  these  equations  mqre  useful  for  real-time  tracking  error 
corrections,  Rowlandson1""  took  advantage  of  the  fact  that  the  error  function 
was  asymptotically  bounded  ls 


From  the  equation 


x  +  (x‘ 


for  the  condition  that 
X2  >>  2 


oo 


1  00  1 

uiw  £  exp  (x3)  f  exp  (-t=>-dt-— +T^+— -,1- 

X 


then  exp  (x2)  [  exp  (-t2)  dt  ~  l/2x 


x 


00 


Since  erf  (x)  =  1  -  j—  J  exp  (-t2) 


dt 


x 


Equations  (9)  and  (12)  with  the  substitution  of  (15)  and  (16)  become 

AR^10'6  nqs.  [1  _ _ &-gg£-U!l 


c  sin  B0  (k  +  g)  exp  [(k  +  g)2] 


] 


and 


tan  B0  (k  +  g)  exp  [(k  +  g)2J 


] 


where  g  =  / cr0  tan  fic 

V  o 


(13) 

(14) 

(15) 

(16) 


(17) 

(18) 
(19) 


and  k  =  Rn  /  c  cos  $ 


(20) 


Since  corrections  to  tracking  errors  must  use  the  measured  radar 
data,  these  equations  were  further  modified  to  replace  R0  by  the  measured 
radar  range,  Re,  hereafter  called  R,  and  the  apparent  elevation  angle,  0O, 
instead  of  the  true  elevation  angle,  80. 


Ray  tracing  data  for  the  NBS-CRPL  surface  corrected  exponential 
atmosphere12  was  used  to  determine  a  new  angle,  Y.  This  new  angle  replaces 
£0  and  was  empirically  developed  to  provide  meaningful  error  data  for  a  wide 
range  of  tracking  conditions,  where 

sin  y  =  sin  0O  +  [K0  +  Ks  exp  (-K3R2)]  exp  (-K10o2)  (21) 
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Except  for  the  propagation  constants,  Ns  and  c,  the  range  and  path 
bending  equations  are  now  a  function  of  the  measured  radar  range,  R,  and 
the  apparent  elevation  angle,  0O,  where 


AR  ~ 


10'6  Ns 

c  sin  y 


10'6  Ns 

tan  y 


[i 


[i 


_ _ g  exp  fg5] _ __i 

(k  +  g)  exp  [(k  +  g)8]  J 

_ g-g£EE_[g!] _ ] 

(k  +  g)  exp  [(k  +  g)?]  J 


(22) 

(23) 


where 


tan  y 


(24) 


and 


cos  y 


y  is  given  by  Equation  (21). 


(25) 


For  station  values  of  Ns  in  the  range  from  340  to  375  N  units,  cor¬ 
responding  to  AFWTR  conditions,  comparisons  with  CRPL  ray  tracing  data 
generated  constants  K0  through  K3,  where 


K0  = 

0. 0175 

(26.  1) 

Ki  = 

245.  0 

(26.  2) 

= 

0.  045 

(26.  3) 

k3  = 

0.  370  x  10"4 

(26.4) 

From  Equation  (14)  it  can  be  shown  that  a  theoretical  restriction  is 
placed  on  the  minimum  value  which  g  can  have,  namely, 


g  = 


tan  y  >  >  J  2 


(27) 


For  elevation  angles  of  four  degrees  and  less.  Equation  (27)  is  far 
from  satisfied,  in  which  case  the  accuracy  of  the  calculations  would  be  in 
error.  However,  careful  selection  of  the  above  constants  K0  through  K3  has 
permitted  the  equations  for  range  error  and  path  bending  to  be  extended  into 
the  small  angle  region. 

Under  long  range  tracking  conditions,  k,  (Equation  (25))  becomes 
much  larger  than  g.  This  situation  causes  a  very  great  simplification  to 
Equations  (21),  (22),  and  (23)  which  then  become 

sin  y  ~  sin  0O  4  K0  exp  (-K^Gq2)  (28) 
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(29) 


AR 


T  w 


IQ'6  Ns 

c  sin  V 

1Q'P  Ns 

tan  y 


(30) 


It  was  shown  by  Rowlandson1’  that  the  elevation  angle  error,  £,  can 
be  calculated  from  the  equation 


R 


£  =  0O  -  sin'1  [~  J  sin  (0O  -  |t|  )  dR] 


(31) 


In  the  long  range  case,  It!,  is  independent  of  range  and  Equation  (31) 
b  ecomes 


€  -  60  -  sin  [sin  (0O  -  r) ] 
or  €  at  T 


(32) 

(33) 


Ray  tracing  results  show  that  c  approaches  r  only  for  ranges  in  excess 
of  3000  km. 


2.  A  Comparison  of  Errors  Calculated  From  Ray  Tracing  Analysis  and 

Rowlandson's  Equations 

For  the  hypothetical  missile  trajectory  postulated  by  Gardner4, 

Figure  20  shows  the  magnitude  of  the  true  elevation  angle,  80,  as  a  function 
of  flight  time  (See  Gardner's  report4  for  additional  tracking  parameters.  ) 
The  radio  refractivity  profile  is  defined  by  a  West  Coast  rawinsonde  profile 
shown  in  Figure  21  It  maybe  observed  that  there  is  a  strong  inversion 
present  and,  therefore,  the  profile  represents  very  non-standard  propagation 
conditions  Superimposed  on  the  profile  is  the  NBS-CRPL  surface-corrected 
exponential  model  profile,  Gardner  uses  a  tracking  station  at  1,  000  feet 
above  mean  sea  level  and  located  on  San  Nicolas  Island,  Therefore,  to  cal¬ 
culate  with  Rowlandson's  equations,  the  station  value  for  Ns  from  the  ex¬ 
ponential  is  340  units.  The  surface  value  for  N  is  found  by  extrapolation  to 
be  approximately  360  units  from  which  the  decay  constant,  c,  is  equal  to 
0.  1548  km'1.12 


During  the  flight  time  from  200  to  475  seconds  the  range  and  tracking 
angles  to  the  vehicle  are  sufficiently  large  to  permit  Equations  (28),  (29), 
and  (30)  to  be  used  to  calculate  the  range  error,  £R,  and  the  total  path  bend¬ 
ing,  T  It  should  be  mentioned  at  this  point  that  all  calculations  and  parameter 
definitions  use  the  metric  system.  However,  for  the  application  to  AFWTR 
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FROM  SAN  NICOLAS 


HEIGHT  (K  FT ) 


FIGURE  21.  WEST  COAST  RADIOSONDE  VERSUS  NBS-CRPL 
EXPONENTIAL 
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tracking  analysis  we  have  generally  presented  the  final  results  in  the  English 
system  of  units. 

Figure  22  shows  a  comparison  of  range  error  calculations  with 
Gardner's  ray  tracing  results.  The  agreement  is  seen  to  be  very  good  down 
to  elevation  angles  of  four  degrees.  At  2.5  degrees  the  error  difference 
corresponds  to  20  feet  at  a  measured  range  of  5.75  x  10'  feet.  Figure  23 
shows  a  comparison  of  elevation  angle  errors,  £.  The  agreement  is  within 
0.  10  milliradians  down  to  four  degrees  elevation  angle.  For  a  burnout  at 
275  seconds,  the  difference  between  curves  represents  a  250  foot  height 
error  at  a  missile  altitude  of  800,000  feet. 

A  second  comparison  is  presented  where  an  East  Coast  rawinsonde 
was  used  to  represent  the  refractivity  profile.  Figure  24  shows  the  rawin¬ 
sonde  profile  with  the  NBS-CRPL  model  profile  superimposed.  In  this  case 
the  tracking  site  was  selected  at  44  feet  above  mean  sea  level  which  cor¬ 
responds  to  Gardner's  analysis  using  the  Point  Mugu  tracking  radar.  This 
analysis  could  equally  well  be  applied  to  the  Eastern  Test  Range  tracking 
system  (MISTRAM)  which  is  located  near  sea  level. 

Figure  25  shows  a  comparison  of  the  range  errors  and  agreement  is 
very  good  down  to  elevation  angles  of  four  degrees.  Figure  26  shows  a  com¬ 
parison  of  the  elevation  angle  errors  with  good  agreement  maintained  to 
elevation  angles  just  under  four  degrees. 

During  the  first  100  seconds  of  flight  the  range  to  the  target  is  too 
small  to  permit  the  use  of  the  simplified  Equations  (28),  (29),  and  (30).  A 
set  of  calculations  over  this  time  period  required  using  the  original  Equa¬ 
tions  (21),  (22),  (23),  and  (31).  Figures  27  and  28  show  a  comparison  of 
calculations  with  Gardner's  ray  tracing  results.  The  West  Coast  rawinsonde 
was  used  to  define  the  propagation  conditions  for  the  ray  tracing  analysis  and 
the  tracking  radar  at  San  Nicolas  Island  (1000  feet  elevation)  was  selected 
for  the  comparative  analysis.  The  elevation  angle  error  differences  at  short 
ranges  and  low  tracking  elevation  angles  is  not  as  good  as  in  the  long  range 
cases.  However,  when  the  range  is  small  one  can  tolerate  a  larger  elevation 
angle  error  since  the  height  position  error  is  the  product  of  the  elevation 
angle  error  and  range 

As  the  tracking  elevation  angle  decreases  the  ray  path  makes  a 
smaller  angle,  9,  with  respect  to  the  inversion  layer.  Referring  to  the  bend¬ 
ing  Equation  (12.  1)  the  magnitude  of  cot  9  and  dN  become  large  as  the  ray 
path  begins  to  graze  the  inversion  layer.  The  increase  in  ray  path  bending 
therefore  becomes  very  large  as  the  ray  passes  through  the  inversion  at  these 
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AND  EXPONENTIAL  PROFILES 
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FIGURE  23.  A  COMPARISON  OF  ELEVATION  ANGLE  ERRORS  FOR 
RADIOSONDE  AND  EXPONENTIAL  PROFILES 


HEIGHT  (K  FT) 


FIGURE  24.  EAST  COAST  RADIOSONDE  VERSUS  NBS-CRPL 
EXPONENTIAL 
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FIGURE  25.  A  COMPARISON  OF  RANGE  ERRORS  FOR  RADIOSONDE 
AND  EXPONENTIAL  PROFILES 
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A  COMPARISON  OF  RANGE  AND  ANGLE  ERRORS  DURING  INITIAL  FLIGHT 


shallow  angles  The  elevation  angle  error,  £,  correspondingly  increases  as 
shown  by  Equation  (33),  Therefore,  in  the  presence  of  large  inversions  such 
as  exhibited  by  the  rawinsonde  profiles  the  elevation  angle  error,  £,  will  in¬ 
crease  very  rapidly  as  the  tracking  elevation  angle  becomes  small  The  ex¬ 
ponential  model  for  refractivity  cannot  represent  refractivity  gradients  of  a 
magnitude  which  commonly  occur  in  an  inversion  and  at  the  same  time 
represent  the  rest  of  the  more-standard  propagation  conditions, 

Gardner’s  summary  and  conclusions  were;  M.  .  A  comparison  between 
actual  extreme  refractive  conditions  for  the  Eastern  Test  Range  (ETR)  and 
the  Pacific  Missile  Range  (PMR)  in  relation  to  ETR  and  NBS  exponential 
atmosphere  was  made  for  a  hypothetical  satellite  launch  as  seen  by  two  widely 
separated  radars  It  has  been  shown  that  for  tracking  elevation  angles  greater 
than  five  degrees,  either  type  of  exponential  atmosphere,  originating  from  a 
measured  surface  index,  can  be  used  for  position  determination  within  100 
feet.  Consequently,  the  only  variable  atmosphere  input  required  for  this 
accuracy  is  the  current  surface  refractive  index  Below  five  degree  eleva¬ 
tion  angle  refraction  corrections  should  only  be  made  with  the  use  of  de¬ 
tailed  rawinsonde  plus  r efr actometer  refractive  index  profiles  using  standard 
ray  tracing  techniques.  11 

Based  on  the  comparisons  of  Rowlandson's  calculations  with  Gardner’s 
ray  tracing  results,  we  concur  that  the  exponential  models  can  be  used  to 
calculate  meaningful  propagation  errors  for  tracking  elevation  angles  down 
to  about  four  degrees  The  position  accuracy  in  range  is  generally  much 
better  than  the  elevation  angle  accuracy  since  the  latter  is  so  greatly  affected 
by  variations  in  the  vertical  refractivity  profile 

The  comparison  also  demonstrates  that  above  five  degrees  tracking 
elevation  angle  the  propagation  errors  can  be  calculated  without  using 
standard  ray  tracing  techniques 

With  reference  to  the  earlier  discussion  on  temporal  effects  it  is 
suggested  that  if  more  detailed  refractivity  profiles  are  required  (for  low 
angle  tracking)  that  rawinsonde  data  be  used  The  rawinsonde  measurement 
would  be  made  in  front  of  the  tracking  radar  and  at  some  convenient  distance 
within  20  nautical  miles 

Ref ractometer  measurements  of  the  refractivity  are  not  recommended 
for  two  reasons.  In  the  first  instance,  the  microwave  r ef ractometer  re¬ 
quires  a  great  deal  of  maintenance  and  calibration  and  without  special  ducting 
facilities  on  the  aircraft  it  cannot  be  used  to  measure  refractivity  in  the 
presence  of  visible  moisture.  Second,  the  time  required  to  make  the  airborne 
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measurements,  then  to  reduce  and  correct  the  data,  would  not  permit  the 
application  of  this  information  in  the  tracking  system  in  near  -  real -time . 

3 §  Comments  on  the  Optimization  of  Rowlandson's  Equations 


A  least  squares  optimization  program  was  used  to  determine  the 
constants  K0  through  K3  to  give  the  best  possible  agreement  with  published 
CRPL.  model  ray  trace  data  1?  Constants  were  determined  corresponding  to 
each  exponential  model  listed  in  the  CRPL  Monograph, 

Figures  29,  30,  and  31  show  the  comparison  for  range  error,  £R, 
total  path  bending,  r,  and  elevation  angle  error,  e,  respectively.  These 
figures  pertain  to  a  surface  refractivity  condition  of  313  N  units.  The  agree¬ 
ment  was  found  to  be  equally  good  for  the  other  model  conditions. 

It  was  also  determined  during  this  analyses  that  the  sin  V  function 
could  be  simplified  without  losing  comparative  accuracy  The  final  form  for 
sin  y  which  will  be  used  hereafter  is 

sin  v  =  sin  0O  +  [K0  +  K2  exp  (-K3R)]  exp  (-K^q)  (34) 

Figure  32  shows  the  values  for  the  optimized  constants  versus  surface  re¬ 
fractivity,  N0. 

The  graphical  results  presented  above  were  derived  with  this  new 
expression  for  sin  y,  replacing  the  earlier  expression  shown  by  Equation  (2l)„ 
In  all  cases  the  differences  between  computed  elevation  angle  and  ray  tracing 
calculations  remain  below  0.  04  milliradians . 

The  excellent  agreement  between  the  sets  of  data  clearly  demonstrate 
that  refraction-induced  radar  propagation  errors  can  be  accurately  calculated 
without  ray  tracing  analysis.  It  is  understood,  of  course,  that  the  equations 
and  constants  developed  in  this  latter  instance  are  designed  to  work  with  the 
CRPL  surface-corrected  exponential  model  atmosphere. 
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FIGURE  29.  A  COMPARISON  OF  RANGE  ERRORS  VERSUS  RANGE 
FOR  RAY-TRACED  AND  ROWLANDSON'S  EQUATION 
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FIGURE  30.  A  COMPARISON  OF  RAY  BENDING  VALUES  VERSUS 
RANGE  FOR  RAY-TRACED  AND  ROWLANDSON'S 
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ROWLANDSON'S  EQUATIONS 
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FIGURE  32.  VARIATION  OF  EMPIRICAL  CONSTANTS  WITH 
SURFACE  REFRACTIVIT Y 


SECTION  V 


AN  ANALYSIS  OF  REFRACTION-INDUCED  ERRORS  USING 
AIRBORNE  REFRACTOMETE.R  J\  EASUREM ENTS  AND  A 
COMPARISON  WITH  THE  CRPL  MODEL 

During  the  period  from  19  January  through  29  May  19t>7,  a  total  of  l6 
airborne  missions  were  flown  in  the  AFWTR  sector  From  these  flights  the 
measurements  made  in  close  proximity  to  the  Tranquillan  Peak  radar  site 
were  analyzed  to  determine  the  magnitudes  of  tracking  errors  produced  by 
refraction  With  reference  to  the  effect  of  spatial  and  temporal  variations,  the 
selected  profiles  nearest  to  the  radar  were  considered  to  give  the  most 
meaningful  description  of  propagation  conditions  affecting  this  radar 

Figures  33  through  36  show  examples  of  the  presentations  used  to 
describe  the  magnitudes  of  elevation  angle  and  range  errors  for  each  measure¬ 
ment  Since  the  total  compendium  of  such  data  is  of  little  general  interest 
the  complete  set  of  such  data  is  not  included  in  this  report. 

In  lieu  of  a  direct  graphical  presentation,  the  errors  were  compiled 
statistically  to  show  the  average  behavior  of  the  elevation  angle  and  range 
errors  together  with  the  standard  deviation  about  the  average 

Figure  37  shows  the  average  elevation  angle  error  behavior  with 
height  for  two  groups  of  test  periods.  The  first  group  covered  the  period 
from  20  January  to  8  February  and  the  second  group  from  28  March  to  29 
May  1967  The  presentation  shows  the  angle  error  data  for  apparent  elevation 
angles  of  100  and  200  milli radians , 

Figure  38  shows  the  average  range  error  variation  with  height  for 
the  two  seasons  and  for  apparent  elevation  angles  of  100  and  200  milliradians 

The  vertical  bars  on  the  curves  define  the  standard  deviation  of  the 
data  about  the  average  values. 

From  reference  19,  the  average  value  of  surface  refractivity  for 
Santa  Maria,  California,  is  listed  as  319.  8  and  322. 7  N  units  for  January  and 
February,  respectively.  Since  Santa  Maria  is  close  to  the  Vandenberg  Base 
and  near  sea  level  (258  feet  msl)  it  is  useful  to  determine  whether  or  not 
these  near-surface  readings  can  be  used  to  generate  a  meaningful  CRPL 
p  rofile 

Using  the  average  Santa  Maria  surface  value  for  January- Feb rua ry  of 
321  25,  the  CRPL  exponential  decay  constant,  c,  is12  0.  1471  km"1.  For  the 
Tranquillan  Peak  radar  at  2500  feet  above  msl  (0.  761  km)  the  CRPL  profile 
gives  a  corresponding  station  value  of  287  N  units. 
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FIGURE  33.  THE  VARIATION  OF  ELEVATION  ANGLE  ERROR  WITH 

HEIGHT  FOR  VARIOUS  APPARENT  ELEVATION  ANGLES 
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FIGURE  34.  THE  VARIATION  OF  ELEVATION  ANGLE  ERROR  WITH 

HEIGHT  FOR  VARIOUS  APPARENT  ELEVATION  ANGLES 
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FIGURE  35.  VARIATION  OF  RANGE  ERROR  WITH  HEIGHT  FOR 
VARIOUS  APPARENT  ELEVATION  ANGLES 


JAN  19,  1967  SPIRAL  A 
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FIGURE  36.  VARIATION  OF  RANGE  ERROR  WITH  HEIGHT  FOR 
VARIOUS  APPARENT  ELEVATION  ANGLES 
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Using  all  available  airborne  r ef ractometer  measurements  made  close 
to  the  Tranquillan  Peak  radar  and  all  Vandenberg  radiosonde  data  for  the 
period  20  January  to  8  February,  the  average  station  value  is  calculated  to 
be  276  N  units.  The  standard  deviation  about  this  average  is  7.  65  N  units. 

It  therefore  appears  that  the  station  value  of  r efractivity ,  276  N  units, 
is  less  on  the  average  then  would  be  determined  from  the  CRPL  exponential 
model  and  a  surface  value  of  about  321  N  units  (Santa  Maria).  Radiosonde 
launch  data  over  the  January-February  period  also  shows  surface  values 
around  320  N  units  from  Vandenberg  and  the  Boathouse. 

L.  Refraction  Errors  Calculated  from  Rowlandson's  Equations 

The  elevation  angle  and  range  errors  were  calculated  using  the  closed 
form  equations  developed  in  Section  IV  (Rowlandson's  equations)  where  sin  v 
is  now  given  by  Equation  (34).  The  appropriate  CRPL  model  is  defined  by  the 
following  method, 

(i)  From  aircraft  and  radiosonde  measurements  the  average 
station  value  is  determined  to  be  276  N  units. 

(ii)  The  index  at  100,  000  feet,  from  standard  tables,  is  constrained 
to  be  3. 8  N  units  for  all  occasions,20 

(iii)  The  simple  exponential  model  becomes  a  straight  line  on  a 
semi-logarithmic  plot,  shown  on  Figure  39.  Projecting  the 
line  from  3  8  N  units  at  100,  000  feet  and  27b  N  units  at  2500 
feet,  an  equivalent  surface  value,  N0,  is  302  N  units. 

(iv)  The  CRPL  exponential  decay  constant,  c,  is  then  0.  140  (km”* 1)12 
for  this  equivalent  surface  value  of  302  N  units. 

The  calculations  are  shown  on  Figures  37  and  38  for  the  elevation 
angle  and  range  errors,  respectively.  The  elevation  angle  calculations  agree 
with  the  average  January-February  results  to  within  0.  1  milliradians .  Re¬ 
ferring  to  Figure  23,  this  agreement  is  essentially  the  same  as  that  obtained 
for  the  hypothetical  missile  trajectory  analysis.  Similarily,  the  range  error 
differences  between  calculations  and  the  January-February  average  (Figure  38) 
is  representative  of  the  differences  obtained  for  the  hypothetical  missile 
trajectory  analysis,  Figure  22.  The  differences  between  the  calculations  and 
ray  tracing  results  become  less  as  the  apparent  elevation  angle,  0O,  is  in¬ 
creased  as  shown  on  Figures  37  and  38 
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This  comparative  analysis  again  shows  that  the  elevation  angle  and 
range  errors  can  be  determined  to  a  high  degree  of  accuracy  with  the  closed 
form  equations 

It  should  also  be  pointed  out  that  the  optimization  of  the  constants  K0 
through  K3  (Equation  (34))  used  the  available  CRPL  ray  tracing  data12  which 
pertains  to  a  radar  at  sea  level.  In  the  above  analysis  the  constants  were 
chosen  for  a  station  value,  Ns,  of  276  N  units  (Reference  Figure  32).  The 
decay  constant,  c,  was  determined  for  a  surface  value,  N0,  of  302  N  units. 

To  obtain  a  final  optimization  of  these  equations  it  will  be  necessary  to  generate 
a  set  of  ray  tracing  data  pertaining  to  a  station  height  of  2500  feet  above  mean 
sea  level  (Tranquillan  Peak  radar)  and  for  propagation  in  a  series  of  different 
CRPL.  model  atmospheres. 

Since  the  analyses  show  very  good  agreement  between  the  data  from 
ray  tracings  and  the  closed  form  equation  we  do  not  feel  that  such  an  additional 
optimization  program  can  be  justified  at  this  time. 
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SECTION  VI 


THE  CORRECTION  OF  TRACKING  ERRORS  AT  VISIBLE 
OPTICAL  FREQUENCIES 


When  visibility  conditions  permit,  it  may  be  required  to  track  vehicles 
with  ballistic  cameras.  Refraction  induced  errors  are  produced  by  variations 
in  the  dry  term  of  the  index  of  refractivity  given  by  Equation  (3). 


NDRY  " 


77.  6  P 
T 


(35) 


It  can  be  shown21  that  in  an  isothermal  atmosphere  the  air  pressure, 

P,  decreases  exponentially  with  height  in  order  to  maintain  hydrostatic 
equilibrium.  A  plot  of  the  variation  of  with  height  would  then  be 

represented  by  a  straight  line  on  the  s emi -logarithmic  plot  of  Figure  40. 

These  conditions  occur  above  40,  000  feet  altitude  as  shown  by  both  the  NACA 
standard  atmosphere2  and  typical  radiosonde  measurements  from  Vandenberg. 


Below  30,  000  feet  the  air  temperature  (degrees  Kelvin)  tends  to  de¬ 
crease  linearly  with  height  in  which  case  the  pressure  decreases  with  height 
according  to21 

T  g/Ry 

P  =  Po  fcr)  <36> 

1  0 

where  p0  =  the  surface  pressure  (mb) 

T0  =  the  surface  temperature  (°  K) 

T  =  the  air  temperature  at  a  given  altitude  (°  K) 
g  =  the  acceleration  of  gravity  (m/ secc) 

R  =  the  specific  gas  constant  (m?/secs/°  K) 

V  =  the  constant  lapse  rate  of  temperature  T  (  K/m) 


Therefore,  the  dry  term  of  refractivity  does  not  follow  a  simple  ex¬ 
ponential  dependence  with  height  over  the  first  40,  000  foot  height  interval. 
Analysis  shows  that  the  dry  term  behavior  with  height  can  be  given  by  an 
empirically  determined  expression 

NDRY  ~  No  exP  [-0.0167  (h1-31)]  (37) 

where  N0  is  the  dry  term  of  refractivity  at  the  surface  and  h  is  the  height. 

In  order  to  use  Rowlandson's  equations  it  would  be  advantageous  if 
the  vertical  variation  of  Nj-^y  could  be  expressed  as  a  simple  exponential 
function  of  height.  A  comparison  was  made  between  the  propagation  errors 
using  the  rawinsonde  measurements  and  a  simple  exponential  model.  The 
exponential  model  is  shown  by  the  straight  line  on  Figure  40  originating  at  a 
surface  value  of  274  N  units  and  passing  through  3. SINT  units  at  100,  000  feet  in 
height. 
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FIGURE  40.  A  COMPARISON  OF  DRY  REFRACTIVITY 
MEASUREMENTS  WITH  MODELS 


Figure  41  shows  the  elevation  angle  error  differences,  Ac,  as  a 
function  of  elevation  angle,  0O.  Ray  tracings  were  carried  out  with  the  target 
at  30,  000  feet  altitude  and  the  tracking  camera  at  100  feet  above  the  surface. 

A  second  comparison  in  Figure  42  places  the  target  at  100,000  feet 
above  sea  level. 


Comparing  Figures  41  and  42,  the  elevation  angle  error  differences 
are  smaller  under  the  longer  range  tracking  conditions.  This  indicates  that 
the  effective  gradients  of  refractivity  which  produce  ray  bending  are  more 
nearly  equal  as  the  height  interval  is  increased.  Referring  to  Figure  40,  it 
is  seen  that  the  magnitude  of  the  radiosonde  gradient  is  less  than  the  simple 
exponential  gradient  certainly  over  the  first  30,  000  foot  height  interval. 
Above  this  height  the  magnitude  of  the  radiosonde  gradient  becomes  larger 
than  the  gradient  of  the  simple  exponential  function.  Therefore,  when  ray 
tracing  over  the  100,  000  foot  height  interval,  the  opposing  radiosonde 
gradients  tend  to  average  out  and  provide  an  overall  effective  gradient  which 
is  closer  to  the  gradient  of  the  simple  exponential  function. 


Figure  43  shows  a  comparison  of  the  range  errors  for  the  radiosonde 
and  the  simple  exponential  profile.  In  this  case  the  total  range  error  dif¬ 
ferences,  A(AR),  increase  as  the  height  interval  is  increased.  The  range 
error,  AR,  is  a  direct  function  of  the  magnitudes  of  refractivity  along  the 
ray  path  where  (reference  Equation  (7)) 


AR 


o 


(38) 


From  Figure  40  it  is  apparent  that  the  absolute  differences  on  the  dry  re¬ 
fractivity  for  the  profiles  in  question  continue  to  be  significant  above  30,  000 
feet  with  the  radiosonde  data  always  larger  than  the  simple  exponential  data. 


Therefore,  the  total  range  error  given  by  Equation  (38)  increases  at 
a  greater  rate  along  the  path  when  radiosonde  data  are  used  than  it  does  when 
the  simple  exponential  is  used.  The  range  error  difference,  A(AR),  there¬ 
fore  increases  directly  as  the  path  length  increases. 

The  magnitudes  of  the  elevation  angle  differences  are  significant  even 
for  tracking  conditions  above  two  degrees  elevation  angle.  Therefore,  the 
use  of  a  simple  exponential  function  to  correct  propagation  errors  in  the 
visible  region  is  not  recommended. 
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FIGURE  43.  RANGE  ERROR  DIFFERENCE  BETWEEN  RADIOSONDE 
AND  EXPONENTIAL  MODEL  (OPTICAL  TRACKING) 


1 .  The  Calculation  of  Optical  Propagation  Errors  With  a  Quadratic 

Exponential  Function 

Freeman  showed  that  a  solution  to  the  ray  tracing  integrals  (Equations 
(7)  and  (12.  1)  could  be  found  by  again  integrating  along  the  slant  path,  R0,  but 
with  a*  quadratic  exponential  function  to  describe  the  vertical  variation  of  N.14 
This  equation  is  given  by 

N  =  Ns  exp  ( -  chs  -  dhs2)  (39) 

where  Ns  is  the  station  value  of  refractivity 
hs  is  the  height  above  the  station  (km) 
and  c  and  d  are  decay  constants. 


Figure  44  shows  the  vertical  variation  of  the  dry  term  of  refractivity 
measured  by  the  radiosonde.  Using  the  quadratic  exponential  function 
(Equation  (39)),  it  is  almost  possible  to  match  the  radiosonde  results.  Care 
was  taken  to  get  good  agreement  near  the  surface  since  the  elevation  angle 
errors  are  very  dependent  upon  the  initial  index  gradient. 

2.  Some  Considerations  on  the  Use  of  a  Bi -Exponential  Atmosphere 

This  section  deals  with  a  refractivity  profile  defined  by  the  sum  of 
two  terms  each  one  expressing  an  exponentially  height-dependent  refractivity 
profile. ^  Let  the  total  refractivity  be  given  by 

N  =  N0  e'Cl  h  +  Nw  e‘C2h  (40) 


The  first  term  represents  the  dry  term  and  the  second  the  wet  term 
of  refractivity.  Figure  45  shows  the  total  value  of  refractivity  obtained  from 
a  radiosonde  (13  June  1967,  SNI).  In  the  discussion  on  optical  tracking  cor¬ 
rections  it  was  shown  that  the  dry  term  of  refractivity  could  not  be  expressed 
as  one  simple  exponential  function  over  the  whole  height  interval  , 


If  the  interval  is  broken  into  two  height  sections,  the  dry  term  can  be 
approximated  by  two  simple  exponentials,  where 


N0  —  280  e 


-0. 0329h 


(0  <  h  <  40  k  ft) 


(41) 


N0- 


-0 . 0508h 

76  e 


(40  ^  h  ^  100  k  ft) 


(42) 


The  wet  term  which  is  plotted  on  the  lower  left  side  of  the  figure  is  approxi¬ 
mated  by  an  exponential 


Nw  —  62 


-0.  1 37  5h 
e 


(43) 
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FIGURE  44.  THE  DRY  TERM  OF  REFRACTIVIT Y 
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FIGURE  45.  A  BI-EXPONENTIAL  REFRACTIVITY  MODEL 


t 


The  contribution  of  the  wet  term  above  30,  000  feet  is  negligible  so 
there  is  no  restriction  on  the  height  interval. 


Using  simple  exponentials  the  total  value  of  refractivity  could  be  ex¬ 
pressed  by 


_ T  ,0n  -0.0329h  -0 .  1 37  5h 

N'j  ~  280  e  +  62  e 


~  76 


•0. 0508h 


(0  £  h  £  40  k  ft)  (44) 

(40  <  h  <  100  k  ft)  (45) 


It  was  shown  previously  that  the  dry  term  over  the  whole  height  in¬ 
terval  could  be  expressed  as 


N  —  Nsd 


-(ch  +  dh2) 
e 


(46) 


Therefore,  another  bi  -  exponential  form  for  the  total  refractivity  could  be 
given  by 


NT-Nsd 


(ch  +  dh2) 


+  N 


-bh 


w< 


(47) 


Using  Freeman's  development  it  is  apparent  that  either  of  the  bi-ex¬ 
ponential  models  can  be  directly  integrated  to  give  an  expression  for  the  re¬ 
tardation  error,  £  R,  and  the  total  ray  path  bending,  T. 14,16  These  equations 
may  then  be  put  into  the  form  of  Rowlandson's  equations  and  empirically 
fitted  to  ray  tracing  data.17 

It  was  already  demonstrated  that  the  simple  exponential  model  for 
total  refractivity  (CRPL  surface  corrected  model)  provides  a  meaningful  de¬ 
scription  of  the  propagation  conditions  for  tracking  elevation  angles  down  to 
around  four  degrees.  Therefore,  at  this  point  there  is  no  justification  to 
further  complicate  the  calculations  through  the  use  of  a  bi  -  exponential  model. 

However,  if  greater  accuracy  requirements  are  found  to  be  necessary 
at  low  tracking  angles  the  bi- exponential  model  permits  greater  flexibility  to 
describe  propagation  conditions  as  measured  by  radiosondes. 
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SECTION  VII 


VELOCITY  MEASUREMENT  ERRORS 


The  apparent  radial  velocity  of  a  target  can  be  determined  by  either 
differentiating  the  range  data  or  by  measuring  the  doppler  frequency  shift  on 
the  returned  carrier  signal.  The  doppler  method  is  generally  preferred  be¬ 
cause'  of  greater  inherent  measurement  accuracy  and  because  differentiation 
tends  to  produce  additional  noise  in  the  measurement. 


Referring  to  Figure  46,  a  target  velocity  V  has  a  component  Vu  in  the 
direction  away  from  the  radar.  Following  the  method  of  Fannin  and  Jehn23, 
let  T  be  the  time  required  for  a  signal  to  travel  from  the  radar  to  the  target 
and  back.  It  can  be  shown34  that  the  doppler  frequency  shift,  £f>  is  related 
to  T  by  the  expression 


dr  _  £f 
dt  ~  f  +  £f 


(48) 


From  Figure  46,  if  u  is  in  the  radial  direction  and  v  in  the  tangential 
direction 


dr 

dt 


*1  .  v  +  — 

6u  u  dv 


Vv 


(49) 


However,  bT/civ  is  zero  since  motion  in  the  tangential  direction  does 
not  change  the  distance  between  the  target  and  the  radar.  If  the  index  of  re¬ 
fraction  is  nt  at  the  target  the  velocity  of  electromagnetic  propagation  is  then 
c/nj-  at  the  target.  Then 


dT 

dt 

2nfr 

c 

'  vu 

(50) 

v  = 

c 

dT 

(51) 

u 

2nt 

dt 

Since  the  true  value  of  n^.  is  generally  unknown,  the  apparent  velocity,  Vua, 
is  then 


V 


ua  2n. 


dr 

dt 


(52) 


The  velocity  error  from  Equations  (51)  and  (52)  is  then 
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Av  =  vua  -  vu 

=  £  <J_  .  J.,  .  51 

2  na  n^  dt 

=  )  •  vu  (53) 

na 


At  high  altitudes  n  is  very  close  to  unity,  to  within  a  few  parts  per 
million.  Therefore,  the  radial  velocity  measurement  error  due  to  the  un¬ 
certainty  in  knowing  the  index  of  refraction  at  the  target  would  be  small  for 
targets  at  very  great  heights. 

The  other  terms  neglected  in  the  above  derivation  is  that  the  round 
trip  time,  T,  will  be  increased  by  the  retardation  error,  AR,  (  reference 
Section  IV),  by  the  difference  in  geometric  path  length,  -  R0,  and  by  the 
fact  that  the  retardation  error,  £R,  is  increasing  in  a  direction  along  the 
ray  path  rather  than  along  R0.  However,  one  can  readily  determine  that 
these  errors  are  negligible  in  comparison  with  the  effect  of  ray  bending 
which  is  now  presented. 

Using  the  methods  of  Millman^6,  it  is  apparent  from  Figure  46  that 
Vr  =  V  cos  (0  +  6)  (54) 

Vu  =  V  cos  if)  (55) 


Due  to  ray  bending  the  target  appears  to  be  moving  along  a  path,  Va, 

where 

v a  =  vu  cos  e  (56) 

The  velocity  error  due  to  ray  bending  is  then 

AV  =  Va  -  Vr  (57) 

Substituting  from  Equations  (54),  (55),  and  (56)  into  (57),  then 

AV  =  V  cos  if)  [—  -  4^-]  +  V  sin  lb  •  sin  6  (58) 

Again,  6s/2  and  e3/2  are  extremely  small,  therefore 

AV  ci  V  sin  ib  *  sin  6  (59) 


The  measurement  of  the  radial  component  of  target  velocity  is,  there¬ 
fore,  principally  in  error  due  to  ray  path  bending  and  is  directly  dependent 
upon  the  tangential  component  of  target  velocity. 
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The  tangential  component  of  velocity  would  be  measured  by  noting  the 
rate  of  change  of  the  true  elevation  angle,  80,  at  particular  ranges,  that  is 

d  en 

V  sin  *  =  R  •  — - ■£- 
dt 

=  R  •  ^  (0O  -  €)  (60) 

Methods  to  compute  e  and  6  have  been  reported  by  Rowlandson.1 

The  doppler  frequency  error  due  to  ray  bending  is  then 
2f 

d  (f^)  =  -  —  *  Y  *  sin  ip  *  sin  6  (6l) 


For  a  spherically  stratified  atmosphere,  and  using  Snell's  Law  in  polar 
coordinates,  Millman  provides  another  method  to  compute  6. 


6  =  cos'1  [^s  COS  (0O  -  O]  -  cos'1  [■ 
rt 


n< 


nt  rt 


cos 


6o] 


(62) 


where  r s 

rt 

6o 

€ 


n 

n 


s 

t 


is  the  radial  distance  from  the  station  to  the  earth's  center 

is  the  radial  distance  from  the  target  to  the  earth's  center 

is  the  apparent  tracking  elevation  angle 

is  the  elevation  angle  error 

is  the  station  value  of  refractive  index 

is  the  value  of  refractive  index  at  the  target. 


For  a  station  height,  Hs,  and  a  target  height,  h,  relative  to  the 
station,  then 

rs  =  r0  +  Hs 

rt  =  r0  +  Hs  +  h 

also  Ns=  1  + 

Nj-  =  1  +  =  1  +  10~6  Ns  e  ^  (for  an  exponential  atmosphere) 

where  Ns  is  the  station  value  of  refractivity 
and  c  is  the  exponential  decay  constant. 

For  most  applications  we  are  concerned  with  herein 


75 


r0  >>h  >>HS 


Ns 

>>“^  >> 
10s 


therefore 

/  ,  h 

rS/rt  “  1  '  — 
r0 

and  r°  ~  1  -  —  +  10'6  (Ns  -  Nt). 

Nt  rt  r0  st 


To  a  first  order,  6,  can  then  be  expressed  by 
6  =  cos'1  [(1  -  h/r0)  cos  (90  -  c)]  -  cos'1  [(1  -  h/r0  +  10'6 
(Ns  -  Nt))  •  cos  60] 


(63) 

(64) 


(65) 


The  measurement  errors  produced  by  refraction  apply  equally  to 
either  doppler  or  range  differentiation  methods  used  to  determine  target 
velocity. 
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SECTION  VIII 


SUMMARY  COMMENTS 

This  report  has  considered  several  aspects  of  the  effect  of  refraction 
on  radar  and  optical  tracking  accuracy.  The  emphasis  was  placed  on  methods 
to  correct  r efraction- indue ed  errors  and  considerable  attention  given  to  the 
use  of  closed  form  equations  for  the  calculation  of  tracking  errors  in  real¬ 
time. 


♦  From  an  operational  point  of  view,  it  was  considered  important  to 
evaluate  the  accuracy  of  these  real-time  correction  methods  because  of  their 
greater  simplicity  over  ray  tracing  methods.  In  particular,  these  closed 
form  expressions  require  a  limited  amount  of  meteorological  support  data 
which  can  be  interpreted  and  used  in  near  -  real- time . 

It  should  also  be  pointed  out  that  ray  tracing  results  are  only  as 
accurate  as  the  meteorological  data  and  unless  exceptional  care  is  taken  in 
the  definition  of  the  propagation  conditions  the  ray  tracing  results  cannot  be 
tacitly  assumed  to  give  an  absolute  measure  of  propagation  errors. 

Below  five  degree  tracking  angles  the  ray  tracing  results  should  pro¬ 
vide  a  more  accurate  description  of  the  propagation  errors  than  the  closed 
form  equations.  This  is  simply  because  the  non- exponential  behavior  of  the 
vertical  variation  of  refractivity  with  height  becomes  more  apparent  at  low 
altitudes.  However,  two  other  considerations  may  impose  a  limit  on  the  use 
of  ray  tracing  methods  below  five  degrees.  In  the  first  case,  any  tendency 
towards  ducting  or  anomalous  propagation  will  begin  to  be  noticed  at  these 
low  elevation  angles.  When  these  effects  are  evident,  ray  tracing  methods, 
as  commonly  defined,  cannot  be  used/*P  ropagation  characteristics  must 
then  be  determined  from  the  field  equations  governing  electromagnetic 
signal  behavior.6  In  the  second  case,  tracking  below  five  degrees  can  sub¬ 
ject  the  radar  to  the  effects  of  multipath  propagation  and  it  then  becomes 
difficult  to  obtain  an  accurate  determination  of  the  apparent  elevation  angle,  0O. 

Most  of  the  ray  tracing  results  used  in  the  comparison  with  the  closed 
form  equations  were  derived  from  radiosonde  measurements.  It  was  shown 
by  Bean  and  Dutton3^  that  unless  the  radiosonde  measurements  are  corrected 

*  for  the  lag  constants  and  the  time  variation  of  air  temperature  on  the  humidity 
sensor,  the  measured  gradients  can  be  significantly  in  error.  Carlson3*7 
shows  that  with  an  uncertainty  in  measuring  the  refractivity  to  within  ±  6  N 
units  in  the  first  1000  feet  and  ±  2  N  units  in  each  higher  1000  foot  layer, 

that  elevation  angle  errors  of  ±  0.  05  milliradians  would  be  experienced  for 
tracking  angles  of  2.  5  degrees.  The  error,  of  course,  falls  off  as  the  eleva¬ 
tion  angle  is  increased. 
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Further  analyses  by  Barnett^8,  Bean  and  Cahoon^,  Mayer30,  Rainey 
and  Thorn31,  Cramond,  et  al32,  Anderson,  et  al33,  and  many  others  clearly 
demonstrate  the  difficulty  in  determining  accurate  and  representative  radio 
refractivity  data  for  use  in  real-time  ray  tracing  analysis. 

To  summarize  these  comments,  I  would  like  to  refer  to  the  paper 
prepared  by  Ratner  and  Bower34  of  the  PMR  Division,  FEC,  VAFB.  Several 
comments  were  made  in  this  report  which  are  pertinent  to  the  analysis  pre¬ 
sented  in  our  (SURC)  report 

(a)  ‘Considerable  investigation  has  been  conducted  at  Point  Mugu 
to  establish  the  experimental  fact  that  the  exponential  atmo¬ 
sphere  model  introduces  sizeable  errors  because  it  is  too 
oversimplified.  11 

(b)  “Radar  tracking  errors  due  to  tropospheric  refraction  can  best 
be  handled  when  the  index  of  refraction  profile  through  the 
troposphere  is  determined.  The  index  of  refraction  profile 
has  a  relatively  systematic  exponential  distribution  with 
altitude,  together  with  a  superimposed  fine  -  structure,  which 

is  not  readily  amenable  to  analytical  treatment.  The  fine- 
structure  is  due  to  local  variations  in  temperature,  humidity 
and  pressure.  The  extent  to  which  tracking  errors  due  to 
refraction  can  be  reduced  is,  in  general,  proportional  to  the 
amount  of  meteorological  data  in  conjunction  with  a  tracking 
operation.  " 

(c)  “REF  RAC  is  used  for  post-flight  data  reduction,  since  its 
execution  cycle  time  is  too  great  to  permit  real-time  analysis.11 

It  has  been  demonstrated  in  this  (SURC)  report  by  a  comparison  with 
Gardner's  data4  that  for  tracking  above  five  degrees  elevation  angle  that  an 
exponential  model  can  be  used  to  determine  refraction-induced  errors  Be¬ 
tween  five  degrees  and  approximately  two  degrees,  directly  measured  radio 
refractivity  data  are  required  together  with  ray  tracing  analysis  to  obtain 
accurate  tracking  error  data.  The  comment,  (a),  above,  is,  therefore,  con¬ 
ditioned  by  the  requirement  to  track  below  five  degrees.  Tracking  below  two 
degrees  is  considered  to  be  impractical  due  to  multipath  propagation  effects 
and  limitations  on  the  applicability  of  conventional  ray  tracing  methods. 

Comment,  (b),  above,  reiterates  the  problem  associated  with  obtain¬ 
ing  accurate  real-time  measurements  of  the  radio  refractivity  structure. 
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This  report  (SURC)  has  attempted  to  show  that  these  fine  -  structure  errors 
can  be  reduced  by  using  real-time  measurements  near  the  site  and  that 
spatial  effects  are  constrained  by  the  effect  of  earth  curvature. 

Finally,  with  respect  to  (c),  the  use  of  Rowlandson's  equations, 
particularly  when  simplifications  of  these  equations  can  be  made  for  long 
range  tracking,  permits  real-time  calculations  of  tracking  errors  to  be 
carried  out  above  five  degree  elevation  angles. 
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SECTION  IX 


RECOMMENDATIONS 

1  A  comparison  of  tracking  error  data  be  made  using  the  AFWTR 
ray  tracing  program  and  Rowlandson's  equations  for  representative  missile 
trajectories,  initially  on  a  post-flight  basis . 

2  The  comparative  data  be  analyzed  with  respect  to  known  tra¬ 
jectory  data,  the  latter  made  available  from  independent  missile  guidance 
and  impact  coordinate  data. 

3.  The  ray  tracing  and  closed  form  (Rowlandson)  methods  to 

derive  tracking  error  data  be  compared  on  the  basis  of  real-time  data  avail¬ 
ability.  This  comparison  will  be  initiated  at  a  pre -designated  time  and 
accurate  records  maintained  of  the  times  required  to  obtain  and  implement 
pertinent  radio  refractivity  data  and  to  process  the  data  for  tracking  correction 
use . 


4.  The  results  of  these  analyses  will  be  presented  to  AFWTR  at 
which  time  a  decision  will  be  made  to  implement  a  technique  or  a  combina¬ 
tion  of  techniques  to  support  real-time  radar  and  ballistic  camera  tracking 
operations . 

5.  The  technique  which  is  selected  will  then  be  re-evaluated 
under  live  launch  conditions. 

6.  Finally,  it  is  recommended  that  additional  studies  be  initiated 
following  the  complete  evaluation  of  these  refraction  correction  methods 
towards  optimizing  the  tracking  data  from  independent  velocity  measuring 
systems . 
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PREFACE 


The  following  paper  describes  a  C-131  Convair  flight  facility  which 
was  instrumented  to  obtain  meteorological  and  radio  refractivity  measure¬ 
ments  . 


The  aircraft,  under  the  direction  of  the  Electronics  Systems  Division, 
USAF,  is  used  to  support  radar  refraction  studies  at  the  National  Ranges  and 
also  meteorological  investigations  for  Air  Force  oriented  programs. 

Having  been  completely  associated  with  this  facility  during  and  since 
its  inception,  my  intention  herein  is  to  provide  a  reasonably  detailed  des¬ 
cription  of  the  aircraft's  capabilities,  its  past  activities,  and  to  indicate  its 
future  potential. 


Lyall  G.  Rowlandson 


A  DESCRIPTION  OF  A  C-  131  CONVAIR  AIRCRAFT 
USED  FOR  ATMOSPHERIC  RESEARCH 


1 . 0  INTRODUCTION 

A  USAF  C-131  Convair  (37812),  based  at  L.  G.  Hanscom  Field, 
Bedford,  Massachusetts,  was  instrumented  in  1962  to  obtain  meteorological 
measurements  associated  with  radio  wave  propagation.  Figure  1  shows  a 

*  photograph  of  the  aircraft  in  flight.  To  determine  propagation  characteristics, 
the  fundamental  parameter  to  be  measured  was  the  radio  refractive  index^^. 
However,  to  obtain  accurate  measurements  it  became  necessary  to  also 
monitor  the  characteristics  of  associated  air  temperature,  pressure,  and 
relative  humidity.  Since  these  measurements  can  be  in  error  due  to  the 
velocity  of  the  air  relative  to  the  sampling  probes,  the  indicated  aircraft 
velocity  and  air  pressure  were  also  required. 

The  requirement  to  collect,  assimilate  and  process  large  quantities 
of  data,  together  with  the  requirement  to  correct  the  free-air  variables  due 
to  aircraft  motion,  led  to  a  digital  recording  system  and  an  associated  data 
analysis  computer  program. 

The  use  of  the  corrected  data  in  propagation  analysis  requires  a 
knowledge  of  where  the  measurements  are  recorded  as  a  function  of  time. 

The  facility  was,  therefore,  further  modified  to  include  an  accurate  radar 
altimeter,  35  mm  telescopic  camera  system  and  a  radio-navigational  position. 

It  is  apparent  that  starting  from  what  initially  appeared  to  be  a  rather 
simple  requirement  led  to  a  much  more  complicated  system.  As  additional 
tasks  were  levied  the  facility  was  expanded  to  include  an  air  sampling  system, 
additional  refractometer s,  tape  recorder,  accelerometers,  drift  meter,  and 
radio  communications  for  long  over-water  operations. 

The  details  of  these  various  instruments  and  programs  will  be  des¬ 
cribed  later.  However,  the  essential  point  is  that  the  facility  has  developed 
into  a  flexible  and  accurate  system  for  investigation  of  parameter  character¬ 
istics  which  affect  radio  signal  propagation.  To  expand  its  mission  one  can 

*  think  of  several  additional  instruments  which  are  available  and  which  would 
permit  greater  involvement  in  atmospheric  measurement  programs;  namely, 
a  weather- radar,  more  extensive  camera  facilities,  three  directional  wind 
velocity  probes  (possibly  acoustic  anemometers),  an  infrared  camera  detec¬ 
tion  system  for  water  temperature  recordings,  air  particle  and  aerosol 
samplers,  lyman-alpha  humidiometer  or  dew-point  hygrometer,  water  drop 
size  spectrometer,  etc. 
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FIGURE  1 


The  primary  power  and  distribution  system  is  adequate  to  meet  these 
further  requirements.  Since  1962,  the  advancements  in  solid-state  develop¬ 
ment  could  lead  to  a  considerable  power,  space  and  weight  saving  in  much  of 
the  existing  equipment.  Certainly  the  aircraft  is  large  enough  to  house  ad¬ 
ditional  equipment. 

1  To  summarize,  the  existing  aircraft  equipment  has  provided  an  ade¬ 

quate  facility  for  lower  altitude  tropospheric  propagation  and  meteorological 
i  investigations.  Its  potential  in  this  important  lower  region  of  the  atmosphere 

could  be  greatly  expanded  by  cooperative  efforts  of  the  various  atmospheric 
science  investigators  to  serve  a  much  broader  field. 


$ 


89 


SPONSORSHIP 


2.  0 


The  C-131  aircraft  is  operated  by  the  United  States  Air  Force  under 
the  direction  of  the  Air  Force  Systems  Command.  It  has  been  based  at  L.  G. 

Hanscom  Field,  Bedford,  Massachusetts,  and  was  assigned  to  support  tech¬ 
nical  investigations  of  the  effect  of  the  troposphere  on  command  and  control 
systems.  Its  experimental  flight  function  was,  therefore,  placed  under  the  * 

Electronics  Systems  Division,  Hanscom  Field,  Bedford.  Historically,  since 
1962,  the  aircraft  has  carried  out  extensive  investigations  at  both  the  National  « 

Ranges  (Eastern  Test  Range,  Cape  Kennedy,  and  Western  Test  Range,  Vanden- 
berg,  California).  Special  joint  programs  have  been  coordinated  with  a  large 
cross-section  of  the  scientific  community  in  those  areas  of  investigation  where 
results  would  have  a  direct  impact  on  the  ESD  mission  and,  at  the  same  time, 
by  agreement,  the  work  could  be  carried  out  without  jeopardizing  the  open 
publication  of  scientific  information . 

These  programs  have  been  extremely  fruitful  in  that  the  Electronics 
Systems  Division  has  gained  valuable  information  and  support  from  the 
scientific  agencies  affiliated  with  the  programs  and,  at  the  same  time,  the 
level  of  competence  of  the  Air  Force-Contractor  personnel  has  been  greatly 
improved. 

As  an  example  of  some  types  of  joint  programs,  I  might  mention  the 
radar-backscatter  experiments  with  AFCRL4,  the  water-vapor  radiometric 
program  with  the  Department  of  Meteorology,  MIT,  the  sea-air  and  anomalous 
propagation  experiments  with  the  Woods  Hole  Oceanographic  Institution6,  the 
multipath-propagation  experiments  with  the  Defense  Research  Board  and  the 
Royal  Canadian  Air  Force6’7*8,  the  DCA  experiments  with  RADC,  ESSA9,  etc. 

The  output  from  these  various  activities  is  reflected  in  the  referenced 
list  of  publications,  not  to  mention  the  associated  papers  which  have  been 
and  are  being  prepared  for  publication  in  the  scientific  journals. 

It  has  been  a  real  credit,  in  my  opinion,  to  the  Electronics  Systems 
Division  of  USAF  and  to  the  administering  project  officers  that  a  reasonable 
degree  of  joint  scientific  effort  has  been  permitted  without  jeopardizing  the 
direct  work  to  the  National  Ranges.  Otherwise,  I  believe  the  program  could  % 

have  become  stodgy  and  limited  and  the  effective  support  to  ESD  in  knowledge 
and  technical  consultation  would  have  been  greatly  reduced  from  what  it  has 
been. 
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DESCRIPTION  OF  THE  AIRCRAFT  INSTRUMENTATION 


The  following  table  provides  a  brief  description  of  the  types  of  in¬ 
struments  which  are  presently  in  the  aircraft,  together  with  their  measure¬ 
ment  accuracy  and  time  constants. 

In  addition  to  the  above  instrumentation,  serious  consideration  is 
being  given  towards  mounting  an  infrared  radiometer  to  measure  sea  surface 
temperature  and  a  weather  radar  to  measure  precipitation  characteristics. 
The  engineering  aspects  of  this  program  have  not  been  discussed  in  any  de¬ 
tail  with  the  Airborne  Engineering  Laboratories  at  Hanscom  Field. 

Figure  2  shows,  from  the  top,  the  AMQ-8  vortex  thermometer,  the 
KS4  aerograph  thermometer  and  relative  humidity  probe,  the  University  of 
Texas  refractometer  cavity  and  the  Rosemount  thermometer  at  the  bottom. 

In  Figure  1,  the  two-cavity,  vertical  N  gradient,  refractometer  maybe  seen 
under  the  starboard  wing  tip. 
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TABLE  I 


Instrument 

Range 

Accuracy 

Time  Constant 

1.  Vortex  Thermometer 

-40°  to  +40°  C 

41/2°  C 

0 .  5  seconds 

2.  KS4  Aerograph 
Thermometer 

-40°  to  +40°  C 

41/ 2°  C 

0 .  5  seconds 

3.  Rosemount  Platinum 

Wire  Thermometer 

-40°  to  +40°  C 

40.  1°  C 

0.01  seconds 

4.  MKS-Pressure 
Transduc  er 

1015  to  500  mb 

±0.  3  mb 

1  second 

5.  Giannini  Air  Speed 
Transducer 

0  to  300  kts 

±0 . 5  kts 

1  second 

6.  KS4- Aerograph 
Relative  Humidity 

0  to  100% 

5% 

1  second 

7.  University  of  Texas 
Ref  rac  tometer  -N 

0  to  400  N 

40.  5  N 

0 .  1  second 
(-3  db) 

8.  Radar  Altimeter 

0  to  20,  000  ft. 

±20  feet 

0,  1  second 

9.  Master  Time-HP 

Clock 

l  in  10s 

0 .  1  second 

(as  normally  set) 

10.  Vertical  Accelera¬ 
tion  -  Strathan 

41. 5  g 

40.  1  g 

0.05  second 

11.  Flight  Research 

35  mm  Camera  on 

A28  Gyro  Mount 

0-20,  000  ft. 

±50  feet 

1  /  1  6  s  econd 

EQUIPMENT  ENGINEERED  BUT  NOT  INSTALLED 


12.  CEC  Dew-Point  (Not  Tested  in  Aircraft  System.  ) 

Hygrometer 

13.  University  of  Texas  Designed  to  measure  N  gradient  over  1  meter 

Ref rac tometer s  -  in  vertical. 

Pod  Mounted  Test  flown  but  further  modifications  required. 
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RECORDING  SYSTEMS  AND  DATA  OUTPUTS 


4.  0 


Three  recording  systems  are  used  on  the  aircraft: 

1  A  36-channel  CEC  galvanometer,  paper  chart,  recorder 

2.  A  l6-channel  Ampex,  FM-analog  tape  recorder 

3.  A  special  purpose  punched-paper  tape  digital  recorder. 

The  prime  data  channels,  1  through  9,  Section  3.  0,  are  recorded  on 
all  systems  The  CEC  chart  recorder  has  a  real-time  developing  unit  which 
permits  the  in-flight  operation  of  these  instruments  to  be  monitored. 

The  magnetic  tape  records  of  the  raw  analog  data  are  extremely 
useful  to  obtain  fine  -  structur  e  and  to  plot  strip  chart  or  X-Y  plots  for  later 
playback  analysis.  Spectrum  analysis  has  been  carried  out  to  a  limited 
degree . 


The  digital  recorder  can  accept  up  to  14  channels  of  information 
every  two  seconds  or  can  be  used  on  decreased  sampling  rates.  The  data 
are  recorded  in  standard  BCD  code,  making  it  very  compatible  for  computer 
pr  oces  sing . 

Figure  3  shows  an  example  of  the  X-Y  magnetic  tape  records  of  the 
variation  of  radio  r ef ractivity ,  N,  with  height.  The  aircraft  is  generally 
moving  in  upward  or  downward  spirals  at  a  500  foot  /minute  rate.  The  hori¬ 
zontal  thickness  on  parts  of  the  record  is  produced  by  the  horizontal  varia¬ 
tions  of  N  as  the  aircraft  orbits  about  the  vertical  axis. 

The  raw  digital  data  obtained  in  punched  paper  tape  are  packed  on 
magnetic  tape  for  direct  input  to  special  purpose  computer  programs. 

Table  II  shows  a  sample  printout  from  the  computer.  The  first  four  sets  of 
data  are  shown,  the  first  group  representing  the  uncorrected  free  air  vari¬ 
ables.  Looking  at  the  second  to  last  column  it  is  apparent  that  data  were  re¬ 
corded  every  two  seconds. 

In  this  particular  test  the  platinum  wire  thermometer  was  not  re¬ 
quired  and  was  not  calibrated  or  used  directly  during  the  flight.  The  aero¬ 
graph  and  vortex  thermometer  were  used  and,  due  to  the  effect  of  air  speed, 
the  aerograph  readings  are  greater  than  the  vortex  thermometer  which  is 
relatively  unaffected  by  speed.  The  relative  humidity  was  high  and  in  such 
a  case  the  carbon  element  in  the  KS4  aerograph  probe  can  indicate  100% 
even  though  the  relative  humidity  may  be  several  percentages  less  The 
radio  r  ef  ractivity ,  N,  together  with  its  variation  from  a  calibration  level 
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ure  3.  .RELATIVE  RADIO  REFRACTIVITY  PROFILE 
DECEMBER,  1965  11:00-13:45  SAN  JUAN  LOCAL  TIME 


SAMPLES  OF  COMPUTER  PRINTOUTS  OF  PROCESSED  AIRCRAFT  DIGITAL  DATA 


INPUT  DATA  FOR  THE  COMPUTATION  OF  FREE  AIR  VARIABLES 

FOR  AIRCRAFT  FLIGHT  -  SAN  JUAN  MISSION  NO.  401,  7  DECEMBER  1965 
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(delta  N)  also  shows  that  the  relative  humidity  was  very  high.  The  last 
column  is  showing  a  manually  induced  event,  number  3.  A  variety  of  event 
numbers  can  be  indicated,  each  of  which  has  some  particular  significance  to 
the  mission.  In  this  case,  event  3  means  llbeginnmg  downward  spiral" 
obviously  from  an  altitude  of  870  meters  above  mean  sea  level 

•  The  second  group  again  shows  some  of  the  initial  free-air  variables 

but  which  have  now  been  corrected  in  the  computer  program  to  take  account 
of  aircraft  speed  and  altitude.  The  agreement  between  the  aerograph  and 
vortex  readings  is  now  much  better.  The  air  pressure  reading  has  also  been 
adjusted  by  1  millibar  from  its  original  uncorrected  value. 

The  computer  next  calculates  the  potential  air  temperature  0,10  the 
dry  term  of  the  index  of  refraction,  N^,  based  on  the  corrected  air  temperature 
and  pressure  readings.  The  fact  that  these  data  are  repeated  twice,  under 
the  heading  aerograph,  indicates  that  these  calculations  were  based  on  the 
aerograph  thermometer  measurements.  Of  course,  any  one  of  the  three 
thermometers  may  be  selected. 

The  next  column  shows  the  mixing  ratio10,  wet  term  of  the  radio  re- 
fractivity,  Nw,  and  the  total,  N,  calculated  from  the  aerograph  thermometer, 
the  aerograph  relative  humidity  probe  (carbon  strip)  and  the  air  pressure 
records.  Using  the  Smith- Weintraub  equation3  for  the  radio  r ef racti vity ,  N, 
the  water  vapor  pressure,  ex0,  and  therefore,  the  mixing  ratio,  r,  can  be 
calculated  since  the  air  temperature,  air  pressure  and  r efractometer -measured 
value  for  N  are  known. 

The  aerograph  calculations  show  a  greater  value  for  the  mixing  ratio 
and  the  N  term  due  to  the  fact  the  carbon  strip  indicated  saturation,  whereas 
the  r efractometer  indicated  below  saturation  conditions. 

The  last  column  shows  the  corrected  total  value  for  N  using  the  re- 
fractometer  but  where  aerograph  thermometer  readings  were  used  to  correct 
for  thermal  expansion  of  the  r efractometer  cavity,  Calibration  of  the  flight 
facility  shows  that  water  vapor  calculations  are  much  more  accurate  using 
the  r efractometer  data  than  with  the  aerograph,  providing  precipitable  water 
is  not  present  within  the  r efractometer  cavity. 

The  final  printout  provides  a  selection  of  parameters  which  have  been 
*  used  to  the  greatest  extent  in  radio  propagation  studies.  The  first  is,  of 

course,  the  corrected  radio  r efractivity,  N,  and  the  geopotential  aircraft 
altitude,  h1  .  A  subsequent  modification  in  the  program  permits  the  radar 
altimeter  measurement  to  be  printed  out  adjacent  to  the  geopotential  height 
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column.  Over  water,  these  radar  measurements  are  extremely  useful  The 
modified  index,  M11,  is  derived  from  N,  the  aircraft  height,  h,  and  an  earth 
radius  factor.  It  is  of  special  interest  in  ray  tracing  analyses11.  The  re¬ 
maining  columns  show  the  potential  air  temperature,  0,  mixing  ratio,  r, 

(using  corrected  r efr actometer  data),  the  air  pressure,  P,  and  the  water 
vapor  pressure,  ei0  (obtained  directly  from  the  Smith- Weintraub  equation) 
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FURTHER  ANALYSES  USING  THE  RADIO-METEOROLOGICAL  DATA 


5.0 
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There  are,  of  course,  many  ways  in  which  the  processed  data  may  be 
used  For  radio  wave  propagation,  the  most  important  information  is  the 
behavior  of  the  radio  r efracti vity ,  N,  in  both  its  vertical  and  horizontal 
dimension 

The  p  receding'  N-  data,  together  with  height,  can  first  of  all  be  used  to 
define  the  vertical  variation  of  N  and  its  effect  on  radar  tracking  errors.  For 
example,  the  following  tables  show  a  comparison  of  errors  using  two  N  pro¬ 
files,  the  first  obtained  with  the  aircraft  and  the  second  obtained  from  a 
radiosonde. 

In  Table  III,  a  typical  set  of  results  are  detei  mined  for  a  radio  ray 
initiated  at  a  point  60.4  feet  above  the  earth’s  surface,  at  an  initial  elevation 
angle  of  18  mr,  and  terminating  at  four  different  stop  heights.  The  profile 
N(h)  was  derived  from  a  particular  aircraft  sounding  in  the  local  area  be¬ 
ginning  at  2018  hours. 

From  the  results  of  the  ray  tracing  analyses  the  electrical  ranges 
are  first  calculated  and  then  the  range  error  produced  by  retardation  of  the 
velocity  of  propagation  along  the  path  (and,  of  course,  due  to  the  N  values) 
relative  to  the  true  geometric  path  See  Figure  4  for  an  illustration  of  the 
parameters  under  consideration 

The  next  column  shows  the  total  amount  of  ray  bending  which  took 
place  and  the  resulting  elevation  angle  error,  E.  The  next  column  shows  a 
k  value  which,  when  multiplied  by  the  true  earth’s  radius,  gives  an  effective 
earth’s  radius,  Ae.  This  effective  radius  represents  an  enlarged  earth  over 
which  one  may  consider  the  radio  waves  to  travel  in  straight  lines11. 

In  Table  IV  the  preceding  results  obtained  with  the  aircraft  measure¬ 
ments  are  compared  with  a  particular  radiosonde  sounding.  The  differences 
may  appear  to  be  small  but  in  many  applications  where  very  accurate  track¬ 
ing  accuracy  is  required  the  results  can  indicate  that  radiosonde  data  are 
not  sufficiently  detailed  to  describe  the  propagation  conditions. 

Similarly,  with  greater  difficulty,  analyses  can  be  carried  out  using 
both  the  vertical  and  horizontal  characteristics  of  radio  r  efracti  vity .  In 
this  case  aircraft  measurements  would  be  required  in  both  the  vertical  and 
horizontal  directions12. 
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TABLE  III 


CALCULATION  OF  REFRACTION  EFFECTS 


^2 

Range 
R  (nm) 

Range 

Error 

(R-Ro)  ft. 

T  otal 

B  ending 

T  (mr) 

Elevation 
Angle 
Error 
€  (mr) 

k 

Effective 

Earth1  s 

Radius 

Ap  (nm) 

8750  0 

58. 135 

98. 84984 

3. 8869643 

1 . 8508379 

1.282458 

4406.  11 

8800. 0 

58  400 

99. 2 1  3'6 1 

3. 9110985 

1 . 8601380 

1  282619 

4406. 69 

8850  0 

58. 665 

99  57561 

3. 9351915 

1 . 8694472 

1 . 282783 

4407  27 

8950  0 

58. 901 

100  29431 

3.9832550 

1 . 8880917 

1 . 283120 

4408  47 

N(h)  =  Profile  I  (A/C  Sounding  Primrose  2018  hours) 
hx  =60.4  feet 

90  =  18  0  mr 


h2  =  Nominal  A/C  Pressure  Altitude,  8850  feet  (Cold  Lake) 
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FIGURE  4.  RAY -TRACING  GEOMETRY 


* 


101 


COMPARISON  OF  REFRACTION  VARIABLES  FOR  TWO  PROFILES 
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In  the  evaluation  of  radiometer  systems,  the  line  integral  of  water 
vapor  and  temperature  are  of  most  direct  interest  but  the  ray  tracing  data 
are  still  required  to  determine  the  path  of  the  signal.  For  examination  of 
anomalous  propagation  conditions  the  vertical  gradient  of  the  radio  re- 
fractivity  is  of  greatest  interest  and  potential  air  temperature  correspond¬ 
ingly  of  greatest  interest  in  atmospheric  stability  analyses11. 

It  is  apparent  that  the  basic  aircraft  data  can  be  used  in  a  wide  range 
of  interconnected  subjects.  Also,  the  addition  of  a  few  pieces  of  selected 
equipment  could  greatly  enhance  the  potential  of  the  facility  in  the  areas  of 
sea-air  interaction,  air  pollution  and  aerosol  characteristics,  rainfall  dis¬ 
tributions  (water  reclamation),  and  weather  modification  programs. 


SUMMARY 
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References  to  the  detailed  data  reduction  and  ray  tracing  programs 
have  not  been  given  since  these  are  of  little  general  interest.  The  essential 
intention  was  to  indicate  the  direct  airborne  measurement  capability  of  this 
aircraft  and  to  suggest  the  potential  it  could  have  for  other  expanded  areas 
of  investigation  in  the  lower  troposphere. 

The  ray  tracing  and  data  reduction  analysis  were  discussed  in  con¬ 
siderable  detail  mainly  because  they  are  interrelated  and  represent  the  im¬ 
portant  output  in  most  radio  propagation  experiments. 

Finally,  the  problem  of  obtaining  accurate  data  with  a  fast-moving 
airborne  platform  is  formidable  and  possibly  was  not  given  sufficient  atten¬ 
tion  This  subject  was  again  considered  to  be  of  little  general  interest  in 
spite  of  the  fact  that  it  is  an  extremely  important  and  difficult  area. 

The  development  of  the  two  sensing  cavity  refractometers  to  measure 
the  vertical  N  gradient  has  been  frustrated  due  to  the  aircraft's  electrical 
noise  generated  on  long  lines  passing  through  the  wing  and  to  microphonics 
associated  with  the  active  RF  elements  However,  recent  measurements 
indicate  that  it  appears  possible  to  measure  0.  1  N  units  variation  between 
the  1  meter  vertically  separated  cavities  New,  open  cavities  are  being 
developed  in  order  to  improve  the  frequency  response  of  the  whole  system13 
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